ADA  124  746 

RADC-TR-82-295.  Vol  II  (of  two) 

Interim  Report 
November  1982 


ACOSS-11  (ACTIVE  CONTROL  OF  SPACE 
STRUCTURES) 

The  Charles  Stark  Draper  Laboratory,  Inc. 


m 


Sponsored  by 

Defense  Advanced  Research  Projects  Agency  (DOD) 
ARPA  Order  No.  3655 


APPROVED  FOR  PUOUC  RELEASE;  DISTRIBUTION  UNLIMITED 


s 


DTIC 

|ELECTE 

FEB  2  2  1983 


The  views  and  conclusions  contained  in  this  document  are  those  of  the 
authors  and  should  not  be  interpreted  as  necessarily  representing 
the  official  policies,  either  expressed  or  implied,  of  the  Defense  Ad¬ 
vanced  Research  Projects  Agency  or  the  U.S.  Government. 


ROME  AIR  DEVELOPMENT  CENTER 
Air  Force  Systems  Command 
Grlffiss  Air  Force  Base,  NY  13441 


*T> 


N 


This  report  has  been  reviewed  by  the  RADC  Public  Affairs  Office  (PA)  and 
is  releasable  to  the  National  Technical  Information  Service  (NTIS) .  At  NTIS 
it  will  be  releasable  to  the  general  public,  including  foreign  nations. 

RADC-TR-82-295 ,  Volume  II  (of  two)  has  been  reviewed  and  is  approved 
for  publication. 


APPROVED 


RICHARD  W.  CARMAN 
Project  Engineer 


*pproved;  (JJljlL 

FRANK  J.  REHM 
Technical  Director 
Surveillance  Division 


Acting  Chief,  Plans  Office 


If  your  address  has  changed  or  if  you  wish  to  be  removed  from  the  RADC 
mailing  list,  or  if  the  addressee  is  no  longer  employed  by  your  organization, 
please  notify  RADC  (  OCSE)  Grif f iss  AFB  NY  13441.  This  will  assist  us  in 
maintaining  a  current  mailing  list. 

Do  not  return  copies  of  this  report  unless  contractual  obligations  or  notices 
on  a  specific  document  requires  that  it  be  returned. 


Volume  II 


ACOSS-11  (ACTIVE  CONTROL  OF  SPACE  STRUCTURES) 

R*  R.  Strunce 
E.  Fogel 
D.  R.  Hegg 
J.  G.  Lin 
J .  D .  Turner 
H.  M.  Chun 


Contractor:  The  Charles  Stark  Draper  Laboratory,  Inc. 

Contract  Number:  F30602-81-C-0180 

Effective  Date  of  Contract:  27  April  1981 

Contract  Expiration  Date:  27  April  1984 

Short  Title  of  Work:  Acoss-11  (Active  Control  of  Space 

Structures) 

Program  Code  Number:  1E20 

Period  of  Work  Covered:  Nov  81  -  Apr  82 


Principal  Investigator: 

Phone : 


Dr.  Keto  Soosaar 
617  258-2575 


Project  Engineer: 

Phone : 


Richard  M.  Carman 
315  330-3148 


Approved  for  public  release;  distribution  unlimited. 


This  research  was  supported  by  the  Defense  Advanced 
Research  Projects  Agency  of  the  Department  of  Defense 
and  was  monitored  by  Richard  M.  Carman  (RADC/OCSE) , 
Griff iss  AFB  NY  13441  under  Contract  F30602-81-C-0180. 


UNCLASSIFIED 


security  classification  or  tHis  page  finiwi  OotomEntotod). 

REPORT  DOCUMENTATION  PAGE 

READ  INSTRUCTIONS 

BEFORE  COMPLETING  FORM 

\.  report1  number  a.  govt  Accession  no. 

RADC-TR-82-295 ,  Vol  II  (of  two) ,4 1) -A 1 

3.  RECIPIENT'S  CATALOG  NUMBER 

4.  TITLE  (end  Subiiil*, ) 

ACOSS-11  (ACTIVE  CONTROL  OF  SPACE 

STRUCTURES)  _ 

S.  TYPE  OP  REPORT  «  PERIOO  COVERED 

Interim  Report 

Nov  81  -  Apr  82 

6.  PERFORMING  O'RG.  REPORT  NUMBER 

CSDL-R-1538  -■ 

7.  Autworfj; 

R«  R.  Strunce  J.  G.  Lin 

E.  Fogel  J.  D.  Turner 

D.  R.  Hegg  H.  M.  Chun 

•.  contract  or  grant  numberoi 

F30602-81-C-0180, 

S.  PERFORMING  ORGANIZATION  NAME  ANO  AOOrESS 

The  Charles  Stark  Draper  Laboratory,  Inc. 

555  Technology  Square 

Cambridge  MA  02139 

10.  PROGRAM  ELEMENT.  PROJECT.  TASK 
AREA  4  WORK  UNIT  NUMBERS 

62301E 

C6550116 

H.  CONTROLLING  OFFICE  NAME  ANO  AOORESS 

Defense  Advanced  Research  Projects  Agency 

1400  Wilson  Blvd 

12.  REPORT  OATt 

November  1982 

t3.  NUMBER  OF  PAGES 

162 

M.  monitoring  agency  NAME  A  AOORESSfff  dllloronl  tram  Controlling  Olllco) 

Rome  Air  Development  Center  (OCSE) 

Griffiss  AFB  NY  13441 

IS.  SECURITY  CLASS,  (of  thie  report) 

UNCLASSIFIED 

ISa.  OECLASSIFICATION/DOWNGRAOING 
N/^cnedule 

IS.  DISTRIBUTION  STATEMENT  (of  this  Roport) 

Approved  for  public  release;  distribution  unlimited. 

17.  DISTRIBUTION  STATEMENT  (ol  tho  abstract  mntmrmd  In  Block  20.  II  different  from  Roport) 

Same 

IS.  supplementary  notes 

RADC  Project  Engineer:  Richard  M.  Carman  (OCSE) 

t9.  KEY  WOROS  (Continue  on  roworoe  «/tf«  if  noceeooey  md  Identify  by  block  ntenber) 

Broad-band  Disturbance  Least  Squares  Approximation 

Power  Spectral  Density  QR  Factorization 

Sensor/Actuator  Selection  Disturbance-Rejection  Control 

i  Node  Shape  Matching  Kalman  Filter 

\  ■  (see  reverse) 

2o^\  ABSTRACT  (Continue  on  reran*  ride  It  no coeoery  ond  Identify  by  block  number) 

Volume  2  of  this  report  describes  progress  in  the  application  of  active 
controller  design  methods  to  the  ACOSS  Model  No.  2  structure.  Attention 
is  focused  on  attenuation  of  the  effects  of  a  broad-band  disturbance 
input  upon  line-of-sight  rotation  error,  and  upon  the  generation  of 
smooth  large-angle  slew  maneuvers.  A  ten-mode  design  model  is  used.  For 
disturbance  attenuation,  a  new  method  for  sensor/actuator  selection  based 
directly  on  line-of-sight  considerations  is  proposed.  Three  disturbance- 

DO  ,  j2S  n  M73  eo'T'o«  or  •  nov  «t  is  obsolete  UNCLASSIFIED 


security  classification  or  this  paoc  f»>w*i  Dim  *»<•*•*) 


UNCLASSIFIED 


fCUWTV  CUAW'ICATIOH  OF  THU  OAOggfcjg  gg»  jggg 


attenuation  designs  are  described.  The  first  assumes  a  complete 
knowledge  of  the  disturbance  statistics,  while  the  others  assume  no  such 
explicit  knowledge.  All  three  designs  indicate  convincingly  that  damp¬ 
ing  augmentation  alone  is  not  sufficient  to  attenuate  effects  of  the 
disturbance;  augmentation  of  stiffness  is  also  required.  Extensions  of 
large-angle  slew  control  to  incorporate  feedback  are  also  presented.  . 

A 

Item  19  (continued)  ' 


Separation  Principle 

Modal  Parameter  Objectives 

Active  Stiffness/Damping  Augmentation 

Augmented- state  Regulator 

Spillover  Reduction 

Actuator  Synthesizers 


Large-Angle  Slew 
Degree-of-Controllability 
Residual  Energy 
Integral  Precompensation 
Riccati  Differential  Equation 


UNCLASSIFIED 


SKCUHITY  CL  AMIAlCATIO*  O'  T“’'  'AOCf»*«»  Of 


ACKNOWLEDGEMENT 


This  report  was  prepared  by  the  Charles  Stark  Draper  Laboratory, 
Inc. ,  under  Contract  F30602-81 -C-01 80.  This  research  was  supported  by 
the  Advanced  Research  Projects  Agency  of  the  Department  of  Defense  and 
monitored  by  the  Rome  Air  Development  Center. 

The  program  manager  is  Dr.  Keto  Soosaar  and  the  project  leader  for 
Active  Control  of  Space  Structures  (ACOSS)  is  Mr.  Robert  R.  Strunce. 

The  technical  director  for  the  ACOSS  work  presented  in  this  report  is 

Dr.  Jiguan  G.  Lin.  The  authors  of  Volume  2  of  this  report  are:  Dr. 

Eliezer  Fogel  (Sections  3,  4),  Dr.  Daniel  R.  Hegg  (Sections  1,  5),  Dr. 

Jiguan  G.  Lin  (Sections  5,  6),  Mr.  Robert  R.  Strunce  (Sections  2,  4), 

and  Dr.  James  D.  Turner  and  Mr.  Ron  M.  Chun  (Section  7).  Assistance 
from  Mr.  Timothy  C.  Henderson  (Sections  2,  3),  Mr.  Warren  J.  Jasper 
(Sections  5,  6),  Mr.  J.  G.  Herther  (Section  5),  Mr.  Robert  R.  Strunce 
(Sections  5,  6),  Mr.  Glen  J.  Kissel  (Section  6),  Dr.  Daniel  R.  Hegg 
(Section  6),  and  Mr.  David  P.  Latimer  (Active  Control  Technology 
Software)  is  gratefully  acknowledged. 


Publication  of  this  report  does  not  constitute  approval  by  the 
Defense  Advanced  Research  Projects  Agency  or  the  United  States  Govern¬ 
ment  of  the  findings  or  conclusions  contained  herein.  It  is  published 
for  the  exchange  and  stimulation  of  ideas. 


i 


Accession  For 

NTIS  GRA&I 
DTIC  TAB 
Unannounced  □ 

Justification _ 


By - 

Distribution/ 


Availability  Codes 


i 

Avail  and/or 

Special 

VOLUME  2 


TABLE  OF  CONTENTS 


Section 


INTRODUCTION  . 

1.1  Scope  . . . . .  i  _  i 

1.1.1  Vibration  Control  .  1_1 

1.1.2  Large-Angle  Slew  Control  .  1-3 

1.2  Limitations  .  !_4 

VCOSS  EXAMPLE  PROBLEM  .  2-1 

2.1  Introduction  .  2-1 

2.2  Broad-Band  Disturbance  .  2-1 

2.3  Input  Disturbance  Implementation  . 2-3 

2.4  Response  to  Input  Disturbance  .  2-5 


SENSOR/ACTUATOR  SELECTION 


3.1  Introduction  . .  3_1 

3.2  The  Node  Shape  Matching  Concept  .  3-2 

3.2.1  Preliminaries  . . .  3_2 

3.2.2  Selection  of  Node  Shapes  .  3-5 

3.2.3  Least  Squares  Formulation  .  3-8 

3.2.4  The  QR  Factorization  . . .  3-9 

3.3  Preliminary  Selection  and  Refinement  .  3-11 


Results 


List  of  References  . . .  3-12 


i  i  i 


Section 


Page 

4  DISTURBANCE  REJECTION  CONTROL  DESIGN  .  4-1 

4.1  Introduction  .  4-1 

4.2  Disturbance  Rejection  Control  .  4-2 

4.2.1  The  Model  .  4-2 

4.2.2  Optimal  Control  Design  .  4-5 

4.3  Application  to  ACOSS  Model  No.  2 . 4-9 

List  of  References  . . .  4-11 

5  LINEAR-QUADRATIC  OPTIMAL  CONTROL  DESIGN  .  5-1 

5.1  The  Design  Problem  . . .  5-1 

5.2  Description  of  Design  Philosophy  .  5-1 

5.2.1  Generation  of  Desired  Closed  Loop  Characteristics  5-2 

5.3  Generation  of  Weighting  Matrices  for  Linear- 

Quadratic  Design  . .  5-9 

5.3.1  Relationships  'Betv/een  Closed-Loop  Modal  Parameters 

and  Weighting  Matrices  .  5-9 

5.3.2  Application  In  Design  .  5-13 

5.4  Design  Results  . .  5-14 

5.5  Numerical  Experience  in  the  Design  Process  .  5-15 

5.5.1  Optimal  vs.  Nonoptimal  Design  . .  5-15 

5.5.2  Conditioning  of  the  Control  Weighting  Matrix  .....  5-15 

5.5.3  State  vs.  Augmented-State  Weighting  .  5-16 

5.5.4  Active  Damping  vs.  Active  Stiffness-Damping 

Augmentation  . . 5-16 

5.6  Work  in  Progress  . . .  5-17 

List  of  References  . .  5-17 

6  LINEAR-QUADRATIC  OPTIMAL  CONTROL  DESIGN  WITH  SPILLOVER 

REDUCTION  .  6-1 

6.1  Introduction  . . . . . . 6-1 

6.2  A  Design  of  Actuator  Synthesizers  . . .  6-1 

6.2.1  Synthesis  for  Prevention  of  Control  Spillover  ....  6-2 


iv 


Section 


6.2.2  Synthesis  for  Correction  of  Modal  Control 

Authority  . . . 

6.3  Optimal  Control  Design  with  Spillover  Reduction 

6.3.1  A  Linear-Quadratic  Optimal  Control  Design  . 

6.3.2  Design  Results  . . . . 

List  of  References  . . 

LARGE-ANGLE  SPACECRAFT  SLEWING  MANEUVERS  . 

7 . 1  Introduction  . . . . . 

7.2  Reduced-Order  Model  Selection  . 

7.3  Retargeting  Maneuvers  For  Flexible  Spacecraft  . 

7.3.1  Optimal  Control  Problem  . 

7.4  Optimal  Integral  Pre -Compensation  LQ  Regulator 


Optimal  Control  Problem  .  7-1 1 

The  Integral  Pre-Compensation  Regulator  Design  ...  7-12 

Some  Frequency  Domain  Features  of  the  Integral 
Pre-Compensation  I£  Regulator  Design  .  7-15 

Observer-Based  Compensation  Design  .  7-15 

Closed-Loop  System  Stability  .  7-17 

A  Closed-Form  Solution  for  the  Differential  Matrix 
Riccati  Equation  for  Fixed  Time  Slewing  Maneuvers  7-19 

Change  of  Variables  for  the  Matrix  Riccati 
Differential  Equation  .  7-20 

List  of  References  . 7-24 


SECTION  1 


INTRODUCTION 


1.1  Scope 

Volume  2  of  the  present  report  gives  an  account  of  the  progress 
made  during  the  reporting  period  in  implementing  various  strategies  for 
active  control  on  ACOSS  Model  No.  2.  The  design  objective  for  the 
several  vibration  control  strategies  employed  is  to  attenuate  the  line- 
of-sight  (LOS)  rotation  error  and  defocus  of  the  optical  system  induced 
by  a  broad-band  disturbance.  A  slewing  controller  for  implementing 
large  changes  in  the  rigid-body  attitude  of  the  structure  with  minimal 
associated  line-of-sight  error  is  also  studied. 

1.1.1  Vibration  Control 

A  general  form  of  disturbance  containing  both  broad-band  and 
discrete-frequency  components  has  been  communicated  by  Riverside 
Research  Institute  to  assist  in  comparing  competing  designs  of  vibration 
controllers.  Several  modifications  of  this  disturbance  specification 
are  made  for  the  purpose  of  the  present  study  in  order  to  reduce 
computation  costs.  First,  the  bandwidth  of  the  broad-band  portion  is  5 
Hz.,  extending  from  0  Hz.  to  5  Hz.  The  reduction  of  the  upper  boundary 
from  15  Hz.  to  5  Hz.  is  dictated  by  the  reduced-order  model  (20  modes) 
being  studied.  The  extension  of  the  lower  boundary  to  0  Hz.  enables  the 
use  of  a  simpler  model  of  the  disturbance  without  sacrifice  of  fidelity 
since  rigid-body  modes  are  not  studied  in  the  disturbance-rejection 
problem.  Second,  the  disturbance  is  assumed  to  be  fixed  in  direction  at 
each  specified  point  of  apj^ication  to  the  structure,  varying 


stochastically  only  in  magnitude.  In  addition,  the  one-sided  spectral 
density  amplitude  in  the  specification  is  retained  for  two-sided  use. 
Details  are  presented  in  Section  2.  These  modifications  do  not  alter 
the  qualitative  objective  of  the  study;  namely,  to  examine  the 
fundamental  difficulties  associated  with  designing  vibration  controllers 
to  stringent  LOS-error  tolerances  in  the  presence  of  a  broad-band 
disturbance. 

A  selection  of  structural  vibration  modes  to  be  considered  in  the 
various  designs,  and  a  corresponding  selection  to  be  used  in  evaluating 
the  designs,  is  determined  by  ranking  each  of  the  156  structural  modes 
in  the  NASTRAN  finite  element  model  with  respect  to  the  root-mean-square 
(RMS)  error  in  LOS-rotation  that  results  from  the  disturbance  acting  on 
that  mode  (viewed  as  a  single  channel).  The  highest  ranking  ten  modes 
according  to  this  criterion  are  chosen  to  constitute  the  model  for  the 
control  designs,  and  the  highest  ranking  twenty  modes  constitute  the 
evaluation  model.  The  upper  limit  on  the  number  of  modes  to  be  included 
is  chosen  with  consideration  of  the  computation  cost.  The  modes 
selected  and  their  ranking  are  also  given  in  Section  2. 

A  choice  of  location  for  sensors  and  actuators  for  the  controller 
design  is  determined  using  an  original  approach  particularly  suited  to 
the  stringent  performance  requirements  in  the  present  design  problem. 
Locations  are  chosen  using  algebraic  methods  which  make  explicit  the 
contributions  to  the  LOS-error  associated  with  particular  sensors  or 
actuators.  Details  are  given  in  Section  3. 

Three  different  approaches  to  the  design  of  a  vibration  controller 
are  presented.  The  first  approach,  disturbance-rejection  control, 
assumes  that  complete  information  on  the  statistical  properties  of  the 
disturbance  is  available.  Using  this  information,  a  stochastic  dynamic 
model  of  the  disturbance  is  concatenated  with  the  open-loop  design  model 
of  the  structure.  A  linear-quadratic-Gaussian  (LQG)  controller  for  the 
enlarged  system  is  designed.  The  design  process  involves  a  few 


iterations  in  which  weighting  matrices  are  adjusted  to  achieve  the 
desired  performance.  Details  are  given  in  Section  4.  Satisfactory 
attenuation  of  the  disturbance  effects  upon  LOS-rotation  error  is 
achieved.  It  is  important  to  note  that  this  attenuation  is  achieved  by 
moving  the  characteristic  frequencies  of  the  closed-loop  system  outside 
the  bandwidth  of  the  disturbance,  as  well  as  by  adding  damping  to  the 
different  modes. 

In  the  second  approach,  linear-quadratic  (LQ)  optimal  control, 
more  realistic  assumptions  with  regard  to  the  knowledge  of  the 
disturbance  characteristics  are  made.  The  disturbance  statistics  are 
not  assumed  to  be  known.  Instead,  the  effect  of  the  disturbance,  as 
reflected  in  the  RMS  values  of  the  LOS-rotation  error,  is  used  to  select 
desired  values  for  closed-loop  modal  characteristic  frequencies  and 
damping  ratios.  Well-known  relationships  connecting  closed-loop  modal 
parameter  values  and  weighting  matrices  for  LQ  design  are  exploited, 
thus  generating  a  systematic  analytical  approach  to  the  selection  of 
weighting  matrices.  Details  are  given  in  Section  5.  As  in  the  first 
design  approach,  results  show  that  LOS-error  reduction  is  obtained 
principally  by  the  augmentation  of  stiffness  rather  than  of  damping. 

The  third  design  approach  is  an  extension  of  the  second  approach; 
actuator  synthesizers  are  incorporated  into  the  design  to  alleviate  the 
problem  of  control  spillover  to  certain  modeled  modes  not  in  the 
controller  design  model.  Details  are  given  in  Section  6. 

1 .1 .2  Large-Angle  Slew  Control 

A  large-angle  slew  maneuver  involves  an  input  of  energy  to  the 
vehicle  which,  in  turn,  induces  vibration  in  the  optical  support 
structure.  The  control  problem  is  to  apply  the  energy  to  the  structure 
in  such  a  way  that  the  desired  attitude  maneuver  is  accomplished  and 
that  the  optical  system  supported  by  the  structure  returns  rapidly  to 
acceptable  levels  of  performance  following  the  maneuver.  Several 
significant  extensions  to  previous  work  are  reported  here.  One  is  to 
solve  the  problem  of  slew  to  a  moving,  rather  than  a  fixed,  target.  The 


other  is  the  development  of  a  closed-loop  control  strategy  using 
integral  compensation  based  on  measurements  which  include  control  and 
control  rate  signals,  rather  than  a  strategy  based  only  on  open-loop 
considerations.  Details  are  given  in  Section  7. 

1 .2  Limitations 

The  present  report  gives  an  account  of  work  in  progress,  rather 
than  of  completed  work.  The  quantitative  results  in  particular  should 
be  viewed  in  that  light. 

With  respect  to  the  vibration  control  designs,  the  present  report 
does  not  include  an  account  of  an  evaluation  of  the  ten-mode  controller 
designs  against  the  twenty-mode  evaluation  model,  or  the  effects  of  the 
use  of  an  observer  in  the  deterministic  LQ  controller  designs.  Investi¬ 
gations  on  these  matters  are  in  progress  and  will  be  reported  on 
subs  eque  nt ly . 

Substantial  work  has  also  been  done  in  other  related  areas,  such 
as  suboptimal  output  feedback,  modal-spring-plus-modal-dashpot  output 
feedback,  and  system  identification,  in  relation  to  flCOSS  Model  No.  2. 

In  addition,  fundamental  investigations  have  been  launched  into  problems 
associated  with  the  explicit  incorporation  of  the  dynamics  of  actuators 
and  of  data-sampling  devices  into  the  mathematical  model  for  the 
structure  and  controller.  Accounts  of  this  work  will  be  reported  in  the 
future. 


SECTION  2 


VCOSS  EXAMPLE  PROBLEM 


2.1  Introduction 

The  problem  posed  for  the  vibration  control  designs  discussed  in 
this  report  is  based  on  the  use  of  Revision  3  (VCOSS  FLEX  Model  No.  2) 
of  ACOSS  Model  No.  2,  more  lightweight  and  flexible  than  the  original 
Model  No.  2.  The  example  problem  is  to  design  a  control  law  for  a 
10-mode  design  model  which  is  driven  by  a  broad-band  disturbance  with  a 
power  spectral  density  (PSD)  of  40  N2*sec  over  a  5  Hz.  bandwidth  and 
to  evaluate  the  resulting  controller  with  a  20-mode  evaluation  model. 
Table  2-1  presents  the  10-mode  and  20-mode  models  and  their  respective 
contribution  to  the  line-of-sight  (LOS)  error  resulting  from  the 
broad-band  disturbance.  The  design  objective  is: 

LOS-X  error  0.05  y-radians 

LOS-Y  error  0.05  y-radians 

LOS-Z  error  _<  25  y-meters 

2. 2  Broad -Band  Disturbance 

The  analytical  two-sided  broad-band  disturbance  is  characterized 
by  a  PSD  which  has  a  second-order  rolloff  (1/s2)  at  5  Hz.  This  is 


modeled  as 


where  G  =  40  N2  «sec 

oi  =  2irf 

c  c 

f  =  5  Hz. 

c 

Eq.  (2-1)  may  be  written  as 


S  (s)  =  F(s )  F(-s ) 

XX 


(2-2) 


where 


s  +  0) 

c 


A  plot  of  the  PSD  for  the  input  disturbance  is  shown  in  Figure  2-1 . 
The  mean  square  value  of  this  disturbance  is 


.2  G  “c  40(2tt5) 

♦  "  "T“  =  - 2 - 


628.32  N2 


This  disturbance  is  consistent  with  the  general  disturbance 
specification  that  has  been  communicated  by  Riverside  Research  Institute 
for  comparison  of  controller  design  strategies. 


2. 3  Input  Disturbance  Implementation 

The  disturbance  sources  are  located  at  node  no.  37  and  node  no.  46 
of  the  VCOSS  FLEX  Model  No.  2.  This  is  true  for  all  modified  versions 
of  the  ACOSS  Model  No.  2  (stiff  and  flexible  models).  The  direction  of 
this  disturbance  makes  an  equal  angle  with  all  three  axes  (x,  y,  z)  at 
each  location.  Figure  2-2  shows  the  geometrical  relationship.  For 
convenience,  it  is  assumed  that  the  input  disturbance  is  fixed  in 
direction  along  the  unit  vector 


U  =  (cos  a)  i  +  (cos  a)  j  +  (cos  a)  k  (2-3) 

where  a  =  54.74° 

The  "projection"  of  the  input  disturbance  on  each  axis  is  the 
disturbance  given  in  Section  2.2. 

With  the  projected  density  Sxx  equal  in  each  direction,  the 
number  of  performance  evaluations  can  be  reduced  by  a  factor  of  three. 
The  actuator  influence  matrix 

node  no.  37  node  no.  46 
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is  multiplied  by  the  matrix 
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to  yield  the  reduced  actuator  influence  matrix 
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Eq.  (2-6)  is  essentially  the  sum  of  the  actuator  influences  in  the  x,  y, 
and  z  directions  for  each  disturbance  location.  This  results  in  only 
having  to  evaluate  two  inputs  with  respect  to  the  line-of -sight 
performance  requirements  instead  of  six  inputs. 


2.4  Response  to  Input  Disturbance 


The  output  response  is 


P  .  .<s) 

out,] 


Hi,i(sl  2  pi„<3)!  i=1'2  •  i-1-2-3 


where  Pin(s)  is  the  PSD  of  the  input  disturbance  described  in  Section 
2.2  and  ,j(s)  is  the  transfer  function  from  the  i-th  component  of 
the  disturbance  location  (node  no.  37  or  no.  46)  to  the  j-th  LOS  error 
function.  Figures  2-3a  through  2-3f  present  the  PSD  of  the  output  for 
the  open-loop  10-mode  design  model.  Figures  2-4a  through  2-4f  present 
the  PSD  of  the  output  for  the  open-loop  20-mode  evaluation  model. 
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Table  2-1 .  VC OSS  example  model , 


RMS  LOS  ERROR 
(micro-radian) 


Sjqc'1'  ■  Input  PSD  In  each  axis 
bjji  ^y»  =  Actuator  influence  in  each  axis 

origin  .  Node  no.  37  or  no.  46 

+  sx*«>  ", +  s«(f)  b*  -  s**m  <wv 


Figure  2-2.  Input  disturbance  geometrical  relationship 
at  node  no.  37  and  node  no.  46. 
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Figure  2-3e.  Disturbance  at  node  46,  output  LOS-Y 
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Figure  2-4b.  Disturbance  at  node  37,  output  LOS-Y 
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SECTION  3 


SENSOR/ACTUATOR  SELECTION 


3.1  Introduction 

The  problem  of  selecting  the  location  and  type  of  sensors  and 
actuators  for  the  active  control  of  large  scale  structures  has  not 
received  enough  attention  by  the  community.  Whereas  such  a  selection  is 
fairly  straightforward  for  well  understood  structures  (e.g.,  beam, 
plate)  it  is  extremely  difficult  when  one  considers  complex  truss 
structures  like  the  ACOSS  Model  No.  2.  To  demonstrate  this  difficulty 
consider  the  utilization  of  member  dampers.  Among  the  advantages  of 
such  actuators  are  their  low  price,  ease  of  installation,  availability, 
and  the  fact  that  they  do  not  affect  the  rigid  body  motion.  However, 
the  number  of  members  in  ACOSS  Model  No.  2  (Revision  3)  is  137,  and  thus 
the  selection  of  location  for  such  actuators  is  formidable.  (Note  that 
member  actuators  can  also  be  placed  between  nodes  which  are  actually  not 
connected  by  members  in  the  original  structure  design,  making  the 
problem  even  more  prohibitive. ) 

In  this  section  we  address  the  selection  of  sensor /actuator 
location  and  derive  a  procedure  which  we  believe  should  provide  a 
meaningful  selection.  Our  approach  to  the  problem  is  demonstrated  via 
application  to  ACOSS  Model  No.  2.  The  selection  provided  here  has  been 
used  in  the  control  designs  discussed  in  the  other  sections  of  this 


A  flow  chart  of  the  selection  procedure  is  given  in  Figure  3-1. 
The  preliminary  selection  step  and  the  "structural  insights"  are 
discussed  in  Section  3.3.  Our  main  contribution,  the  LOS  matching 
concept,  is  the  subject  of  Section  3.2. 


3. 2  The  Node  Shape  Matching  Concept 
3.2.1  Preliminaries 

To  establish  convenient  notations  and  terminology,  a  few 
definitions  are  called  for  at  this  point. 

The  vibration  model  of  a  structure  can  be  specified  via  the 
triplet  (S,  Z,  ♦)  where 

S  =  [a>1 ,  u>2,  •••]  denotes  the  structural  frequencies 

Z  =  ;2,  •••]  denotes  the  associated  damping  ratios 

♦  =  f  »  •  •  •  1  denotes  the  matrix  of  mode  shapes 

-1  -2  J 

(Note:  *  is  the  modal  matrix  associated  with  Z  =  0. )  Each  of  the 

vectors  J.  is  a  6N  vector  which  describes  the  i-th  mode  shape  at  the 

N  structural  nodes  in  6  degrees  of  freedom.  Denote  the  rows  of  $  by 

£2'  i,e*' 


* 


T 


(3 


where  ^  denotes  the  transpose  of  ♦.  Assume  without  loss  of 
generality  that  the  f^  are  organized  such  that  if 


(i) 


4-1 


mod  6 


j-1,  <i-j+l>/6 


(3 


c: 


then  f^  corresponds  the  4-th  mode  in  the  j-th  degree  of  freedom. 

By  analogy  to  the  concept  of  mode  shape,  we  define  node  shape  as 
follows: 

Definition.  Let  x  be  any  vector  of  dimension  6N.  Then  is  called 
a  node  shape  associated  with  x. 

To  illustrate  the  usefulness  of  the  node  shape  concept,  some 
examples  are  given: 

(a)  The  Standard  Node  Shape 


The  vector  f^  corresponds  to  the  node  shape  at  the  4-th  node  in 
the  j-th  degree -of -freedom  using  the  relations  (3-2).  Thus  choosing  the 
vector  x  to  have  the  entries  i=l,  •••,  N,  satisfying 

1a  ,  i  =  64  +  1 

0  ,  i  =  64  +  2 

y  ,  i  =  64  +  3 

0  ,  otherwise 


with  a2  +  02  +  y2  =  1  gives  the  normalized  translation  node  shape 


in  the  direction  0  at  node  4. 

_Y_ 

(b)  LOS  Node  Shape  (See  Ref.  3-1  for  details) 

The  line-of-sight  (LOS)  error  about  the  x  and  y  axes  is  a  linear 
combination  of  the  displacement  of  the  mirror  support-points  expressed 


LOS  =  l  a.  u,  , 
i=1  1 


where  NS  =  number  of  mirror  support  degrees-of-f reedom 


Since  tu  =  n  »  the  above  can  be  expressed  as 


NS 


LOS  =  l  a.$.  n 
i=1 


r -T  yr 
L®  aJ  J1 


where 


— 

T 

“losx, 1 

“losx, NS 

^LOSX 

T 

“losy,i 

“losy,ns 

2losy 

T 

“defocus,! 

“defocus, NS 

_ 2defocus_ 

to  giVe  -LOSX  =  *  ^OSX 


-LOSY  *  ?LOSY  '  -DEFOCUS  “  *  ^DEFOCUS 

the  three  LOS  errors. 

(c)  Axial  Node  Shape 

The  axial  displacement  of  a  member  connecting  nodes  1  and  2  is 


given  by 


where 


A^  ^  =  axial  displacement  between  nodes  1  and  2, 

Uj  =  the  displacement  vectors  at  nodes  1  and  2, 

respectively,  of  the  connecting  element  (1,2) 
b1 ,  b2  =  coordinate  transformation  for  element  (1,2). 

Expanding  this  equation  to  include  all  elements  of  interest,  we  get 
A=*BU,  where  B  is  the  matrix  of  element  coordinate  transforms.  Transfor¬ 
ming  the  displacement  vector  into  modal  coordinates  gives 


A 


T  T  T 
(4  B1) 


n 


and  thus  4TBT  gives  the  node  shape  matrix  for  axial  displacements. 

Each  column  of  the  matrix  represents  the  axial  displacement  of  one 
element  for  a  unit  amplitude  of  each  mode. 

In  light  of  the  above  definition  and  examples  we  note  that 
selecting  the  location  and  type  of  actuators  (e.g.  member  damper, 
momentum  wheel)  and  sensors  is  equivalent  to  choosing  a  set  of  node 
shapes  in  the  sense  of  the  above  definition.  We  propose  the  following 
philosophy  for  the  selection  of  a  set  of  node  shapes. 


3.2.2  Selection  of  Node  Shapes 

The  actuators /sensors  location  selection  can  be  formulated  using 
the  following  dynamic  representation  of  the  structure 


n  +  seh.  +  n  .1 


(3-3) 


P  =  (4  Xp)  n 


v  =  (4‘xv)‘  n 


(3-4) 


where 


SB  =  diag  (2?  u  ,  25  01.,  •••  lx,  u>  ) 
1  1  z  z  n  n 

2  2  2 

52  =  diag  (uj,  ,  u)_,  •••  ui  ) 

12  n 


u.  is  the  input  (actuation)  vector  of  dimension  tt^ 
jd  is  a  position  output  (sensing)  vector  of  dimension  nip 
_v  is  a  velocity  output  (sensing)  vector  of  dimension 
d  is  a  disturbance  input  to  the  structure. 


The  matrices  X  (n  x  m  ),  X  (m  x  n),  and  X  (m  x  n)  generate  the 
u  u  p  p  v  v 

node  shapes  associated  with  the  corresponding  actuators/sensors 

location.  For  example,  if  x^O )  is  the  first  column  of  Xu,  then 

T 

4  xu<1)  is  the  node  shape  associated  with  u-| — the  first  actuator 

in  the  vector  u^.  Similarly,  the  node  shape  associated  with  o-j ,  the 

T 

first  position  sensor,  is  4  x  (1),  where  x  (1)  is  the  first  column 

-p  -p 


of  Xp. 


Noting  that  LOS  =  we  can  formulate  the  selection  problem 


as  follows: 


Find  Xu,  Xp,  Xv  to  enable  the  best  control  of  the  LOS. 
Namely,  if  a  certain  mode  does  not  affect  the  LOS  significantly 
(corresponding  to  a  relatively  small  entry  in  f* )  one  does  not 
have  to  control  it  or  observe  it. 


This  idea  can  be  expressed  mathematically  in  view  o£  Bqs.  (3-1)  and 
(3-2)  as  follows. 

Select  Xp  (and  Xy  and  Xy)  such  that 

•  T 

f  =  4  X  •  a 

-*  P  - 

where  oi  is  some  arbitrary  vector  of  appropriate  dimension  and  =  means 
"approximately  equal". 

The  above  formulation  can  be  somewhat  simplified,  assuming  an 
initial  (large)  set  of  locations  has  already  been  selected  (see  next 
section  for  this  step).  Let  the  corresponding  node  shapes  be 
{f.|,  •••#  f.j}.  The  problem  thus  reduces  to  a  selection  of  a  subset 

of  L  node  shapes  out  of  the  N  which  provides  a  good  approximation  of 
f*. 

Criteria  for  selection 

Choose  nodes  {i-j,  i2*  •••  ijJ  from  {l,  »•»  n}  such  that 

L 

( 1  )  7  o.f .  *  f*  ,  and 

•  ,  3-i. 

j=1  3 

(2)  All  the  oj's  are  of  comparable/acceptable  order  of  magnitude. 

Definition  A  selection  of  node  shapes  satisfying  the  criteria  (1)  and 
(2)  is  called  a  balanced  selection. 
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Remarks 


(i)  The  objective  is  to  achieve  a  balanced  selection  with  L  as  small 
as  possible.  A  balanced  set  of  sensors  (actuators)  means  that 
similar  sensors,  i.e.,  with  the  same  sensitivity,  can  be  used  at 
all  locations  (the  same  power  is  required  of  all  the  actuators). 

(ii)  It  is  obvious  that  a  set  of  independent  nodes  are  sought,  barring 
redundancy  (fault  tolerance)  considerations. 

We  propose  as  a  mathematical  criterion  for  "approximately  equal" 
the  quadratic  cost  function.  Thus  the  selection  problem  can  be 
formulated  as  follows. 


3.2.3  Least  Squares  Formulation 
Solve 


NT  N 

mjn  -  l  ai  D  [**  -  l  «i  <3-5> 

-  i=1  i=1  1  1 

where  D  =  diag  [dj,  d2»  •••  djj],  subject  to  the  constraints: 

(a)  N-L  of  the  a^'s  are  zero,  and 

(b)  The  L  nonzero  a^'s  form  a  balanced  set. 

To  the  best  of  our  knowledge,  the  above  problem  has  not  been 
solved.  Thus,  we  propose  an  iterative  procedure  to  obtain  an  acceptable 
(suboptimal)  approximation  to  the  selection  problem  which  is  based  on 
the  £R  factorization. 

The  purpose  of  the  D  matrix  in  Eq.  (3-5)  is  to  weight  the 
different  modes  according  to  their  relative  contribution  to  the  LOS 
performance.  For  example,  in  ACOSS  Model  No.  2,  with  the  disturbance 
discussed  in  Section  2,  a  reasonable  way  to  weight  the  modes  is  in 
proportion  to  the  root-mean-square  (RMS)  LOS  amplitude  associated  with 
the  given  disturbance. 


Consolidating  the  set  of  node  shapes  ,  *•*, 


into  a  matrix  F, 


and  each  x  denotes  an  arbitrary  value.  Using  the  QR  factorization,  we 
also  have  the  additional  property 


rii  i  ri.t,  i+i  i  • 


and  each  r^  is  larger  than  any  other  entry  in  its  row.  Using  the  QR 
factorization,  Eq.  (3-6)  can  be  written  as 


R  a 


L.S. 


R  a 
u- 


h  i 


(3-8 


«  =  Ru 


Note  that  ^(c^f^)  is  the  projection  of  f#  on  and  that  r^q^  is  the 
length  of  the  component  of  f^,  orthogonal  to  the  span  of  {f^,  •••,  f^ 

Thus  the  pair  (r^ ,  b^)  gives  an  indication  of  the  contribution  of 

f .  to  the  L.S.  approximation  of  f  .  Thus,  elimination  of  f. 's  with 
•1  “1 


relatively  small  r^  or  b^  or  r^*b^  will  not  result  in  a  significant 


increase  in  the  error  of  the  L.S.  approximation,  i.e.,  in  the  norm  of 
the  residual  vector  R 


3.3  Preliminary  Selection  and  Refinement 

As  discussed  in  the  previous  section,  the  selection  procedure  is 

simplified  by  starting  with  an  initial  set  of  locations,  say 

f  »  •••.  f  • 

-1 '  '  -N 

The  preliminary  N  actuators/sensors  have  been  selected  as  member 
actuators  (sensors).  This  selection  is  based  on  the  ranking  of  the 
members  according  to  their  elongation  for  each  mode.  The  members  under¬ 
going  the  most  elongation  (stress)  for  those  inodes  which  have  most 
effect  on  the  LOS  are  selected  as  the  initial  set. 

Application  of  the  selection  procedure  of  the  previous  section 
leads  to  a  reduced  set  of  members  with  LOS  node  shape  matching  which 
does  not  differ  considerably  from  the  match  obtained  with  the 
preliminary  selection.  If  the  matching  is  not  satisfactory,  additional 
actuators  have  to  be  added.  Inspection  of  the  mismatch  in  the  residuals 
as  defined  in  Eq.  (3-9)  indicates  the  modes  which  are  not  well-matched 
(i.e. ,  controlled  or  observed  from  the  selected  set.)  The  actuators 
(sensors)  to  be  added  should  compensate  for  this  deficiency.  As 
discussed  in  the  next  section,  it  is  possible  that  there  exists  no 
member  capable  of  controlling  certain  modes  which  do  have  significant 
influence  on  the  IOS.  In  such  an  event,  a  different  type  of  actuator 
(force  or  torque  generating  actuators)  must  be  utilized.  The  selection 
of  the  locations  for  such  actuators  can  be  aided  by  visual  inspection  of 
the  vibrating  structure  simulation  and/or  the  analysis  of  the  mode 
shape. 


3. 4  Results 

The  selection  procedure  described  above  has  been  applied  to  ACOSS 
Model  No.  2. 

Test  1 .  A  50  mode  model  has  been  used.  The  initial  set  of  actuators 
included  all  mirror  mounts  as  well  as  some  members.  The  selection 
procedure  has  been  applied  to  yield  a  perfect  match  of  the  LOS  where  the 


locations  selected  were  the  mirror  mounts.  This  academic  test  verifies 
the  basic  philosophy. 

Test  2.  For  the  purpose  of  the  control  design  demonstration,  the 
20-mode  model  has  been  selected.  The  weighting  matrix  D  of  Eq.  (3-5) 
has  been  selected  so  that  each  mode  is  weighted  by  its  expected  RMS 
error  due  to  the  prescribed  disturbance.  An  initial  set  of  member 
actuators  has  been  selected  using  the  elongation  concept.  Iterating  the 
selection  procedure  resulted  in  unsatisfactory  residuals.  Examination 
of  the  problematic  modes  revealed  that  the  uncontrolled  modes  were  solar 
panel  modes  which  have  strong  coupling  to  the  LOS.  Moreover,  it  was 
observed  that  there  exists  no  member  capable  of  providing  good  control/ 
observation  authority  on  these  modes.  Adding  force  actuators  proved  to 
be  satisfactory.  The  final  selection  is  given  in  Figure  3-2.  Note  that 
the  selection  is  not  obvious  intuitively,  but  can  (after  the  fact)  be 
given  acceptable  structural  justification.  Moreover,  it  is  conceivable 
that  by  changing  the  properties  of  the  selected  members,  the  structure 
will  be  less  sensitive  (in  terms  of  LOS  performance)  to  the  given 
disturbances. 
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SECTION  4 


DISTURBANCE  REJECTION  CONTROL  DESIGN 

4.1  Introduction 

The  broad-band  disturbance  described  in  Section  2  poses  a 
difficult  control  problem.  Specifically,  the  RMS  value  of  the  line-of- 
sight  (LOS)  due  to  the  disturbance  is,  in  the  open  loop,  of  the  order  of 
100  p-radians  (prad)  in  LOS-X  and  over  1000  prad  in  LOS-Y  (see  Table 
4-1),  whereas  the  desired  RMS  of  the  LOS  error  should  not  exceed  0.05 
prad.  Namely,  the  control  objective  is  to  reduce  the  effect  of  the 
disturbance  by  four  orders  of  magnitude!!! 

To  assess  our  ability  to  achieve  this  design  objective,  an  optimal 
disturbance  rejection  controller  is  presented  in  this  section.  The  goal 
here  is  to  assess  the  following  issues:  (a)  Can  the  control  objective 
be  achieved  under  optimal  (unrealistic)  conditions?  (b)  If  the  answer 
to  (a)  is  affirmative,  what  insight  into  the  design  of  a  realistic 
control  might  be  gained  from  the  characteristics  of  the  optimal  control? 

It  is  demonstrated  in  this  section  that  the  control  objective  can 
be  achieved  under  optimal  conditions,  and  conclusions  are  drawn  with 
regard  to  the  type  of  control  that  should  be  designed  to  reject  the 
given  level  of  disturbance. 

In  Section  4.2,  the  control  design  philosophy  is  outlined. 

Section  4.3  is  devoted  to  a  detailed  description  of  the  designed 
controller.  The  frequency  response  of  the  controlled  structure,  the 
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feedback  gains,  and  the  pole  distribution  of  the  closed  loop  system  are 
presented.  We  also  make  an  effort  to  extrapolate  from  this  design  to 
the  desired  features  of  a  realistic  design. 

4.2  Disturbance  Rejection  Control 
4.2.1  The  Model 

The  theory  of  optimal  disturbance  rejection  is  most  conveniently 
studied  via  the  state-space  representation  of  the  system  dynamics 
(structure  and  disturbance  model)  [Ref.  4-1],  Thus,  we  use  the 
following  representation  for  the  N-mode  structure 

•  2N  m  P 

x  =  Ax  +  Bu  +  B,d,  xeR;ueR,deR  (4- 

u  d 

o 

y  =  Cx  ,  y  e  R  (4- 

LOS  =  C_x  ,  LOS  e  R3  (4- 

Ii 


where,  as  is  common  in  the  context  of  structures,  A  is  given  by: 


®u  =  actuators  influence  matrix  of  dimension  2N  x  i, 
m  =  number  of  actuators 

C  *  sensors  observation-matrix  of  dimension  t  x  2N, 

£  =  number  of  sensors. 

Here  we  consider  colocated  sensors /actuators,  thus  £  =  m.  The 
disturbance  influence-matrix  is  B^,  where  in  our  analysis  it  is  of 
dimension  2N  x  2.  (See  Section  2. ) 

Finally,  since  the  objective  is  to  control  the  structure  subject 
to  the  disturbance,  as  specified  via  its  power  spectral  density  (PSD)  in 
Section  2,  it  is  convenient  to  consider  the  disturbance  as  a  stochastic 
process  represented  by  the  state  equation 

d  =  A  d  +  B  a)  (4-5) 

d  til 

where 


Here,  Xi  and  X2  reflect  the  bandwidth  of  the  disturbances  and  b^, 
b^  their  amplitudes. 

Combining  the  structure  representation  of  Eqs.  (4-1 )-(4-4)  with 
the  disturbance  model  of  Eqs.  (4-5) -(4-6)  yields  the  augmented  system 
equations 


or,  equivalently,  with  proper  definitions  of  the  augmented  terms  in  Eq 
(4-7),  the  augmented  system  dynamic  equations  can  be  written  as 


x  =  A  x  +  B  oj  +  B  u  (4-8a) 

a  a  a  aw  au 


y 


C  x 
a  a 


( 4-8b ) 


LOS 


'aL 


x 

a 


(4-8c) 


The  model  of  the  system  in  Eq,  (4-8)  is  the  basis  for  the  control 
design  below.  It  is  thus  appropriate  at  this  junction  to  question  the 
availability  of  such  a  model  for  a  control  design.  Specifically,  one 
should  question  the  validity  of  the  disturbance  model.  The  assumption 
that  the  spectral  content  of  the  disturbance  is  available  is  a 
reasonable  one,  as  the  disturbance  source  is  often  understood  or  can  be 
measured.  Hence,  assuming  knowledge  of  the  terms  defined  in  Eq.  (4-6) 
is  reasonable.  The  difficult  assumption  is  that  made  with  regard  to  the 
knowledge  of  the  matrix  Bd.  Although  certain  noise  sources  (e.g. , 
imbalance  of  a  rotating  wheel)  have  well-defined  points  of  contact  with 
the  structure  (i.e.,  well-defined  Bd),  other  noise  sources  may  have 
more  distributed  effect  on  the  structure,  such  that  knowledge  of  Bd  is 


less  likely.  The  principal  value  of  the  control  design  study  to  be 
described  below  is  in  the  exposition  of  the  characteristics  and 
limitations  of  a  disturbance  rejection  controller  which  is  based  on  the 
assumed  knowledge  with  regard  to  the  disturbance. 


4.2.2  Optimal  Control  Design 

Given  the  augmented  dynamic  equation  Eq.  (4-8),  we  design  the 
optimal  disturbance  rejection  controller  using  the  separation  principle 
[Ref.  4-1]: 

(a)  We  assume  that  the  augmented  state  is  available,  and  design 
an  optimal  regulator. 

(b)  Using  the  available  measurement  y,  we  construct  a  Kalman  type 
filter  to  estimate  the  state.  It  should  be  noted  that,  in 
our  problem  formulation,  there  is  no  measurement  noise.  We 
use  the  noise  variance  as  a  design  parameter  which  helps  in 
the  loop  shaping  and  circumvents  the  problem  of  designing  a 
"singular  Kalman  filter"  or  am  arbitrary  observer. 

(c)  We  cascade  the  state  estimator  by  the  optimal  state  feedback 
gain. 

The  design  of  the  optimal  regulator  involves  the  following 
procedure. 

(a.1)  Define  the  cost  functional 


T 

^  /  [ (LOS)TQ(LOS)  +  p  UTRU]dt 


J 


lim 


(4-9) 


Free  design  parameters  in  Eg.  (4-9)  are  the  matrix 


and  the  scalar  p. 

We  have  arbitrarily  fixed  R  to  be  the  identity  matrix,  since  we 
assume  that  all  actuators  are  of  the  same  type;  i.e.,  their  power 
consumption  and  influence  are  the  same,  such  that  they  should  be 
penalized  equally.  Obviously,  if  some  actuators  are  more  expensive  to 
use  in  some  sense  (e.g. ,  a  propellant  actuator),  this  should  be 
reflected  in  the  corresponding  entry  of  the  matrix  R.  The  cost 
functional  depends  on  the  line-of -sight  output  variable  "LOS"  and  not 
explicitly  on  the  state,  since  we  are  interested  only  in  line-of-sight 
performance.  The  free  parameters  q^,  i=1,2,3,  are  chosen  to  reflect 
the  contribution  of  the  i-th  degree  of  freedom  to  the  line-of-sight 
error.  The  scalar  p  is  an  adjustable  parameter  to  be  tuned  to  achieve 
the  design  objectives. 

(a. 2)  The  regulator  gain  G  is  calculated  via  the  equations: 


-1  T 

G  =  -  (pR)  BP 

cl 


AT  P  +  PA  +  (C  T  )T  Q(C  _  )  +  PB 

A  a  al,  al,  a 


,(pR) 


-1 


T 

B  P 


( 4-1 Oa ) 
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(4-1  Ob) 


The  regulated  system  (with  state  feedback)  is  given  by 


x  *  (A  +  B  G)  x  +  B  oj  (4-1  la) 

a  a  au  a  aa) 


LOS  =  C  Tx 
aL  a 


(4-1 1b) 


(b)  The  optimal  state  estimator  is  calculated  using  the  relations  dual 
to  those  used  in  the  optimal  regulator  design.  The  estimator  equation 
is 


d 


x 

a 


A  x  +  B  u  +  K  [y  -  C 
a  a  au  1  a 


where  the  optimal  estimator  gain  K  is  calculated  via 

T  -1 

K  -  L  C  (pi) 

a  O 

T  T  T  -1 

AL+LA+B  B  -LC(pI)  CL=0 
a  a  au  au  a  o  a 


(4-12) 


(4-1 3a) 


(4-1 3b) 


The  gain  calculated  in  Eq.  (4-13)  is  optimal  with  respect  to  the  cost 
functional 
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T  * 
x  )  ( x 

a  a 


(4-14) 


where  E  denotes  the  expectation  operator.  The  scalar  pQ  is  a  design 
parameter  in  this  section. 


0 


In  optimal  estimation  theory,  pQ I  is  replaced  by  the  measurement 
noise  covariance  matrix. 

Note;  The  Riccati  equation  Eq.  (4-1  Ob)  is  of  order  2N+p.  Since  the 
disturbance  state  d  is  uncontrollable,  the  Riccati  equation  can  be 
partitioned,  to  yield 


*  P11A  *  5  -  P1,BuR‘'BuP11  -  0  I4-'5a) 


(A  -  B  R_1B  P.jV.  +  P, _A  +  P..B.  =  0  (4-1 5b) 

u  u  11  12  12  d  lid 


where  Eq.  (4-1 5a)  is  a  Riccati  equation  of  order  2N  (quadratic  equation 
in  Pi  i  )  and  Eq.  (4-15b)  is  known  as  a  Sylvester  equation  [Ref.  4-2]. 

The  optimal  gain  is  given  by 


K  -  -  (p*)-’  (b*  P„  B*  P12] 


( 4-1 5c ) 


For  a  large  number  of  disturbances  (p  large;  see  Eq.  (4-1)),  solving  Eq. 
(4-15)  rather  than  Eq.  (4-10)  has  definite  numerical  advantages. 

However,  in  our  case,  p  «  2N  which  reduces  the  benefit  in  solving  Eq. 
(4-15),  especially  when  we  note  the  additional  software  needed  to  solve 
the  Sylvester  equation.  Nevertheless,  it  is  important  to  note  the 
relations  in  Eq.  (4-15),  since  it  clearly  shows  that  the  stability  of 
the  state  feedback  does  not  depend  on  the  disturbance  model.  The  latter 
enters  into  the  plant  only  in  the  feed-forward  term 
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Although  partitioning  of  Eq.  (4-13)  does  not  yield  the  same 
benefit  as  in  the  optimal  state  feedback,  it  should  be  noted  that  the 
stability  of  the  disturbance  rejection  scheme  still  does  not  depend  on 
the  accuracy  of  the  disturbance  model.  This  follows  from  the  separation 
theorem  [Ref.  4-1 ]  which  implies  separate  pole  assignment  for  the 
regulator  and  the  estimator.  Namely,  N  of  the  poles  are  assigned  by  the 
matrix  K,  as  if  the  filter  were  not  in  the  loop,  and  N  poles  are  the 
filter  poles.  Thus,  from  the  stability  point  of  view,  the  questions 
raised  above  concerning  the  disturbance  model  are  of  no  consequence. 
Inaccuracy  of  the  disturbance  models,  as  discussed  above,  will  only 
affect  the  RMS  value. 


4.3  Application  to  ACOSS  Model  No.  2 

The  disturbance  rejection  control  scheme  outlined  above  has  been 
applied  to  ACOSS  Model  No.  2.  The  design  involved  few  iterations  in 
which  the  various  design  parameters  mentioned  earlier  have  been  adjusted 
and  the  corresponding  performance  (RMS  values  and  frequency  response) 
assessed.  It  was  found  that  the  control  penalty  has  to  be  of  the  order 
10“17  to  achieve  the  design  objectives.  Namely,  to  achieve  four 
orders  of  magnitude  reduction  in  the  RMS  value  requires  substantial 
control  effort  which  translates  into  lew  control  penalty  in  the 
quadratic  cost  function  Eq.  (4-9).  Furthermore,  best  results  have  been 
obtained  when  LOS-Y  was  weighted  10  times  more  than  LOS-X,  whereas  LOS-Z 
(defocus)  can  be  lightly  weighted  in  the  cost  penalty.  A  sample  of  the 
RMS  LOS-error  for  different  design  parameters  is  given  in  Table  4-1 . 

For  each  design,  values  are  given  both  for  the  case  that  state  feedback 
is  applied  directly  to  the  state  (LQR)  and  to  the  estimated  state 
( OBS/LQR ) .  In  Figures  4-1  a  through  4-1 f,  the  closed-loop  frequency 
response  with  optimal  state  feedback,  and  the  open  loop  response  from 
u^,  i=l,2,  (ui  =  disturbance  at  node  37;  U2  =  disturbance  at  node 
46)  to  yj,  j=1 ,2,3,  (LOS-X,  LOS-Y,  LOS-Z)  are  plotted.  Figures  4-2a 
through  4-2f  give  the  corresponding  closed-loop  responses  with  observer 
in  the  loop,  and  give  a  comparison  with  the  responses  for  the  closed- 
loop  sy  tem  without  an  observer.  All  figures  are  for  Case  5  of  Table 


Examination  of  Table  4-1  and  Figures  4-1  a  through  4-2f  reveals  the 
following  features  of  the  disturbance  rejection  controller. 

There  is  a  reduction  of  the  gain  over  the  0-5  Hz  range  of  150  db 
to  200  db.  There  is  some  difficulty  in  controlling  the  three  closely- 
spaced  inodes  13,  14,  15  (0.59  Hz,  0.61  Hz,  0.63  Hz  respectively),  a 
phenomenon  which  exhibits  itself  in  the  "glitch"  at  about  0.6  Hz  in 
Figures  4-1a,  4-1c,  4-1d,  4-1e  and  4-1f.  (See  explanation  in  Ref.  4-3.) 
However  this  is  a  local  phenomenon  which  does  not  affect  the  RMS  value 
significantly.  Figures  4-2a  through  4-2f  are  the  corresponding 
frequency  responses  when  a  Kalman  filter  is  introduced  into  the  loop. 
There  is  good  correspondence  of  the  frequency  response  in  each  case 
between  the  optimal  closed  loop  with  and  without  the  filter  in  the  low 
frequencies  (less  than  5  Hz).  For  higher  frequencies  the  optimal 
control  with  filter  has  higher  gain  to  yield  a  broader  bandwidth 
transfer  function.  Note  in  particular  the  deterioration  in  performance 
for  the  channel  u-j  to  y2  where  the  transfer  function  is  almost  flat 
with  the  filter  in  place,  resulting  in  increase  of  the  RMS  LOS  value 
from  0.002  urad  to  0.130  yrad.  A  possible  explanation  of  this  severe 
deterioration  in  performance  is  the  difficulty  in  estimating  the  state 
of  the  disturbance  at  node  37.  This  difficulty  might  stem  from  the 
sensors  location,  since  only  one  sensor  is  located  on  the  upper  support 
truss  compared  to  six  on  the  lower  support  truss  and  two  on  the  solar 
panels  (see  Section  3).  Node  37  is  on  the  upper  support  truss.  Thus, 
changing  the  sensor  configuration  might  rectify  this  difficulty. 

Further  insight  to  the  behavior  of  the  optimal  control  can  be 
gained  by  examining  the  shift  in  the  pole  positions  due  to  the  feedback 
(Table  4-2).  It  is  observed  that  although  some  damping  is  added  to  the 
various  modes,  significant  effort  is  made  in  shifting  the  modes  to 
higher  frequencies,  i.e.,  out  of  the  disturbance  range.  This 
observation  is  supported  also  by  observing  that  the  regulator  gains 
(Table  4-3)  are  of  the  order  10®  for  position  feedback  (columns  1-10) 
and  only  of  the  order  10^  for  the  velocity  feedback  (column  11-20). 

Note  that  relatively  low  gains  are  applied  to  the  disturbance  states. 
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We  conclude  from  the  above  observations  that  achieving  the  design 
objectives  under  realistic  conditions  requires  sophisticated  control 
design.  Specifically,  it  was  observed  that  damping  augmentation  will 
not  suffice  even  if  a  high  level  of  damping  were  possible.  The 
difficulty  in  controlling  closely-spaced  modes  also  indicates  the 
possible  collapse  of  control  schemes  (ours  or  others)  if  some  weakly 
excited  modes  are  ignored  in  the  proximity  of  controlled  modes. 
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Table  4-2.  Effect  of  disturbance  rejection  control  on  poles  distribution 
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Figures  4-1 .  Output  power  spectral  density  from  closed- 
loop  linear  quadratic  regulator  (LQR) 
(LOS-output  formulation). 
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Figures  4-2.  Comparison  of  output  power  spectral  densities 

of  closed-loop  linear  quadratic  gaussian  regulator 
(OBS/L£ft)  with  linear  quadratic  regulator  (LQR). 
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Figure  4-2c.  Disturbance  at  node  37,  output  LOS-Z. 
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Figure  4-2f.  Disturbance  at  node  46,  output  LOS- 


SECTION  5 


LINEAR -QUADRATIC  OPTIMAL  CONTROL  DESIGN 


5. 1  The  Design  Problem 

Preceding  sections  have  described  a  flexible  structure  model,  and 
a  general  model  of  a  disturbance  containing  both  broad-band  and 
discrete  frequency  components.  A  reduced-order  subsystem  (design  model) 
is  defined  which  consists  of  the  ten  modes  that  display  the  largest 
mean-square  line-of -sight  (LOS )  rotation  error  when  the  disturbance  is 
applied  to  the  full  open-loop  structural  model.  Using  the  same 
criterion  for  ranking,  a  twenty-mode  subsystem  for  evaluation  of  the 
design  (evaluation  model)  is  also  defined.  The  problem  treated  in  the 
present  section  is  to  design  a  linear-quadratic  regulator  for  the 
ten-mode  reduced-order  subsystem  which  is  asymptotically  stable  with 
respect  to  the  larger  twenty-mode  reduced  order  subsystem,  and,  in 
addition,  which  assures  that  the  steady-state  LOS  rotation  error  is 
within  the  specified  tolerance  in  spite  of  the  disturbance. 


5.2  Description  of  Design  Philosophy 

The  disturbance-rejection  controller  presented  in  Section  4 
assumes  full  knowledge  of  the  statistics  of  the  disturbance.  With  such 
an  assumption,  it  is  possible  to  include  a  dynamic  model  of  the 
disturbance  as  part  of  the  model  of  the  reduced-order  system  to  be 
controlled.  In  contrast,  for  the  purpose  of  the  present  procedure,  no 
information  about  the  disturbance  is  used  explicitly  in  the  design 
process.  Instead,  the  LOS  errors  generated  by  the  disturbance  acting  on 
the  open-loop  system  are  translated  into  desired  values  for  modal 
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damping  ratio  and  characteristic  frequency  of  the  closed-loop  ten-mode 
subsystem.  If  such  values  cam  be  attained,  a  closed-loop  system  that 
meets  the  overall  design  objectives  in  stability  and  performance,  at 
least  with  respect  to  itself,  if  not  with  respect  to  the  evaluation 
model,  can  be  expected.  These  desired  parameter  values,  which 
indirectly  reflect  certain  effects  of  the  disturbance,  are  used  as 
explicit  objectives  of  a  linear-quadratic  design.  Well-known 
connections  between  such  parameter  objectives  and  the  weighting  matrices 
in  the  quadratic  cost  functional  are  exploited. 

The  approach  just  outlined  embodies  a  systematic  analytical 
procedure  for  selection  of  weighting  matrices  that  is  explicitly  and 
directly  connected  with  the  system  performance  specifications.  It 
avoids  to  a  substantial  degree  the  trial-and-error  process  often 
associated  with  the  selection  of  weighting  matrices  for  linear-quadratic 
regulator  design. 

5.2.1  Generation  of  Desired  Closed-Loop  Characteristics 

First  the  transfer  matrix  connecting  the  LOS-error  outputs  and 
the  disturbance  inputs  for  the  open- loop  structure  is  generated.  Next, 
using  spectral  density  concepts,  the  mean-square- value  of  these  errors 
due  to  the  disturbance,  as  a  function  of  open-loop  modal  parameters,  is 
developed.  Finally,  these  relations  are  used  to  generate  implicit 
relations  for  desired  closed-loop  modal  parameters  by  imposing  an 
allowable  tolerance  on  LOS-errors. 

5. 2. 1.1  Transfer  Matrix:  Disturbance  Inputs  to  LOS-Error  Outputs 

The  open- loop  structural  model  in  modal  coordinates  with 
disturbance  inputs  is: 

n  +  2ZQn  +  n2n  *  (*tb  (5-i) 
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LOS 


where  n  =  (ni » •  •  •  >nij)T  is  the  modal  coordinate  vector,  N  being  the 
number  of  modes  to  be  controlled,  ft  =  diag  (wj , . . .  ,W|j):  NxN  is  the 
matrix  of  modal  characteristic  frequencies,  Z  =  diag  (Ci, . . . , £N): 

NxN  is  the  matrix  of  modal  damping  ratios,  *  =  [♦1|  •••!♦**]*  6vxN  is 
the  matrix  whose  columns  are  the  mode  shape  vectors,  Bq:  6vx$i  is 
the 

structural  influence  matrix  of  the  disturbances,  u„  =  ILw,,:  6,  xl 

D  d  D  1  „ 

T 

is  the  disturbance  input  to  the  structure,  yT __  =  (y  ,,...,yTne  „) 
is  the  vector  of  LOS-error  outputs,  and  C^gg:  £x6v  is  the  combination 
matrix  defining  the  LOS-error  outputs.  The  actuator  equivalence  matrix 
Rg:  6i xS2  reflects  simplifying  assumptions  on  the  disturbance 
allowing  a  lower-dimensional  disturbance  vector  wg:  62x1  (cf. 

Section  2.3;  esp.  Bq.  (2-5)).  Superscript  "T"  denotes  matrix  transpose 

Since  the  system  of  Eq.  (5-1)  is  diagonal  in  the  modal  coordi¬ 
nates,  and  Bq.  (5-2)  can  be  written  in  the  form: 
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it  makes  sense  to  consider  each  single-mode  channel  connecting  the 
disturbance  input  wD  to  the  channel  output  y^QS*  The  corresponding 
transfer  matrix  Ha(s):  £x52  is 
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5. 2. 1.2  Mean-Square-Value  of  LOS -Error  Outputs 

It  is  assumed  that  the  components  of  the  input  disturbance  vector 
wjj  are  uncorrelated  wide-sense-stationary  random  processes.  The 
associated  spectral  density  matrix  is  thus: 

S  (to)  =  S  (to)  I. 
w„  I  on 


where  Sj(to)  is  a  second  order  approximation  to  the  specified 
disturbance  density: 

A  Gto^ 

S  (to  )  -  _ c  ,  -00  <  to  <  +  00 

2  2 
(0  +  to 
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and  Ig  is  the  &2X&2  identity  matrix.  Values  for  the  amplitude 
G  and  bandwidth  frequency  toc  are  specified  in  Section  2.2. 

The  mean-square-value  matrix  for  the  IOS-error  output  random 
process  is  obtained  by  integrating  the  output  spectral  density  over  all 
frequencies  [Ref.  5-1): 


E  a  a  (t) 

yL0S  '  yL0S 


1  1 

=  —  f  S  a  (to)  dto  =  — —  /  S  a  (-is)ds 

9  ir  t  v  2iri  '  v _ 


21,1  ii»  yL0S 


(5-3) 


where 


Sot  (to)  =  H(i<o)S  (to)H(-ito). 
yL0S  °  “D  ° 
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Eq.  (5-3)  reduces  to  the  evaluation  of  a  single  scalar  complex  integral: 


/  F  (s  )F  ( -s  )ds ) 


or  2iri 
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where 


A  2  T 
M  =  G  a)  K  K 
a  c  a  a 


F(s)  = 


(a)  +  s)L  (s) 

c  a 


The  complex  integral  in  Eq.  (5-4)  may  be  evaluated  using  the 
residue  theorem;  convenient  tabulations  of  such  integrals  with  F  a 
general  rational  function  have  been  compiled  (Ref.  5-2],  The  result,  as 
a  function  of  the  modal  parameters,  is: 
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The  root-mean-square  (RMS)  value  for  the  errors  associated  with 
LOS  rotation  is  of  principal  concern  and  is  defined  by: 
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where  the  are  the  upper-left  diagonal  elements  of  the  mean-square- 
value  matrix  in  Eq.  (5-4).  It  follows  that 
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where  ma  is  a  constant  independent  of  the  modal  frequency  and  damping 
ratio. 
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5. 2. 1.3  Implicit  Relations  for  Desired  Closed-Loop  Modal  Parameters 

The  RMS  LOS-rotation  error  induced  by  the  input  disturbance  acting 
on  the  a-th  mode  open- loop  subsystem  is: 


aO 


=  [ RMS  1 

L  LOS,aJ0 


(5-7) 
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where  oJqQ  is  the  characteristic  frequency  and  SgQ  is  the  assumed 
modal  damping  ratio.  For  each  of  the  ten  modes  selected  to  be 
controlled,  the  values  in  Eq.  (5-7)  exceed  the  allowable  value  of  0.05 
g-radians.  A  desired  value  [RMSLos,a]d  for  t*ie  ^  LOS  rotation 
error  is  chosen,  and  we  seek  a  strategy  for  compensation  of  the 
structural  open-loop  system  leading  to  corresponding  parameters 

of  the  closed-loop  compensated  system  which  achieves  the  desired 
error  reduction,  i.e.. 


p  .  =  [ RMS  ].  =  ifu  .,C_,)1/2  ®  (5-8) 

ad  1  L0S,aJd  ^  oed  ok3;  a 

Equation  (5-8)  does  not  reflect  any  modal  coupling  effects  that  result 
from  compensation.  Combining  Eqs.  (5-7)  and  (5-8)  leads  to  an  implicit 
relation  for  the  parameters  ua(j  and  5^.  In  view  of  the  form  of 
Eq.  (5-5),  it  is  convenient  to  express  the  implicit  relation  in  terms  of 
new  variables  y,  k  defined  by: 
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which  has  two  real  roots  of  opposite  sign.  Using  the  positive  root, 

(Yd)*,  one  obtains  the  desired  closed-loop  modal  parameter  pair 

ft  a)  .)  =  [(Y.J./Za,  a)  sought.  Choosing  a  =  u>  corresponds  to 
'ad  ad  'd+  '  au 

leaving  the  a-th  mode  natural  frequency  unchanged. 
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Special  Case  2: 


<  Y*  <  +- 


b  >  0 


This  corresponds  to  fixing  the  closed-loop  damping  ratio  at 
the  value  1/2b,  without  further  constraining  the  natural  frequency 
Equation  (5-9)  becomes  a  quintic  equation  in 


+  (2Vob')Yd  +  (2Vob')yd  f  ('^Yd  +  ^‘“c^  =  0  (5_11) 
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which  has  at  least  one  positive  real  root.  If  (yd)  +  is  a  positive 

real  root  of  Eq.  (5-11),  then  the  closed-loop  system  with  modal 

parameter  pairs  ft  ,,w  ,]  =  fl/2b,  b(Y„)  1  satisfies  the  RMS  LOS- 
v  ad  ad"  N  d  +' 

error  specifications.  Choosing  b  =  1/25aQ  corresponds  to  leaving  the 
a-th  mode  damping  ratio  unchanged. 

Results  of  calculations  performed  to  determine  modal  parameter 
objectives  for  the  design  are  shown  in  Table  5-1 .  The  first  set  of 
three  columns  shows  the  open-loop  modal  parameters  and  the  RMS 
LOS-rotation  error  that  result  from  applying  the  disturbance  to  the  full 
open-loop  structural  model.  The  second  set  of  three  columns  shows  the 
closed-loop  modal  damping  ratios  that  would  be  required  to  reduce  the 
RMS  error  in  LOS  rotation  to  the  allowable  level  of  pad  =  0.05 
p-radians  without  changing  the  modal  frequencies.  These  correspond  to 
solutions  of  Eq.  (5-10)  with  a  =  <aQ q.  Implementation  of  these  ratios 
is  not  desirable,  since  each  represents  real  poles,  one  of  which  would 
have  a  very  large  time  constant  (i.e.,  damping  ratio  far  in  excess  of 
critical).  The  next  set  of  three  columns  shows  that  a  much  more 
reasonable  set  of  modal  damping  ratios  is  obtained  by  raising  the 
allowable  RMS  error  to  pad  =5.0  y-radians.  The  parameters  listed  in 
the  final  set  of  three  columns  reflect  changes  associated  with  the  first 
four  (overdamped)  modes,  whose  parameters  are  adjusted  so  as  to  fix  the 
desired  damping  ratio  at  the  "optimal"  value  of  1//2.  These  correspond 
to  solutions  of  Eq.  (5-11)  with  b  =  1//2.  These  parameter  values  are 
used  in  subsequent  calculations. 
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It  should  be  observed  that  the  implicit  relation  of  Eq.  (5-9) 
allows  considerably  more  flexibility  in  design  than  is  actually 
implemented  with  Special  Cases  1  and  2.  These  special  cases  impose 
restrictions  on  variations  in  the  free  parameters  (pacj*  Ya »  ic^) 
for  the  sake  of  convenience,  as  indicated.  However,  it  is  quite 
conceivable  that  there  may  be  an  advantage  in  dropping  such  restrictions 
and  focusing  directly  on  pole  assignment  instead  of  attempting  explicit 
assignment  of  characteristic  frequencies  and  damping  ratios. 

5. 3  Generation  of  Weighting  Matrices  for  Linear-Quadratic  Design 

In  the  preceding  subsection,  an  approach  to  generating  desired 
values  for  the  closed-loop  poles  was  developed  so  as  to  ensure  that  the 
effects  of  the  disturbance  on  the  LOS  rotation  are  kept  at  an  acceptable 
level.  In  the  present  section,  these  parameter  objectives  are  used  to 
guide  the  selection  of  weighting  matrices  for  a  full-state-feedback 
linear-quadratic  controller  design. 

5.3.1  Relationships  between  Closed-Loop  Modal  Parameters  and  Weighting 
Matrices 

The  open-loop  structural  model  in  modal  coordinates  with  actuator 
inputs  is: 


"ri  +  2Zfln  +  ft  n  = 
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where  BA:  6vxm  is  the  structural  influence  matrix  of  the  actuators, 
and  u  =  (ui,...,um)T  is  the  actuator  input  to  the  structure. 
(Remaining  notation  is  as  in  Eqs.  (5-1)  and  (5-2)).  Defining  state 
variables  by 
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leads  to  a  set  of  N  systems  of  state  equations  of  the  form 


s  =  A  s  +  B  w 
a  a  a  a  a 
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For  each  a,  consider  the  problem  of  minimizing  the  functional 
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along  the  trajectories  of  the  system  Eq.  (5-13),  with  "full-state 
feedback  of  the  form 


F  s 
a  a 
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Note  particularly  that  the  effects  of  coupling  between  inodes  is  not 
modeled  by  the  functionals  in  Eq.  (5-14). 

The  connection  between  the  modal  parameters  of  the  second-order 

closed-loop  system  corresponding  to  Eqs.  (5-13)  and  (5-15)  and  the 

weighting  matrices  in  Eq.  (5-14)  is  well-known  (Refs.  5-3,  5-4].  These 

connections  have  been  exploited  in  multivariable  applications  [Ref. 

5-5].  Essentially,  one  simply  solves  the  Riccati  equation  associated 

with  the  optimization  problem  analytically  and  makes  the  appropriate 

c  c 

connections.  Denoting  the  closed-loop  modal  parameters  by  ui  ,  £  , 

a  a 

the  results  are: 
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where 


Since  the  off-diagonal  elements  of  do  not  enter  into  Eqs.  (5-16)  or 
(5-17),  there  is  no  loss  of  generality  (at  the  level  of  the  second-order 
system  (5-13))  in  assuming  that  is  diagonal. 

An  appropriate  quadratic  cost  functional  for  the  full  N-mode 
design  model  that  is  consistent  with  Eqs.  (5-16)  and  (5-17)  is  the 
following  augmented-state  regulator  formulation: 


J  -  2«Va  4  r  7*5,  t  yLOS  *  ■  *  nw  ♦  »*■]* 


where  the  augmented-state  vector  y  ia  defined  by: 


ya  '  yLOS,«  '  “■  ' .  * 
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The  augmented-state  weighting  matrix  Q  in  Eq.  (5-18)  has  the  structure: 


Ql  0  0 

e  =  o  qd  o 

0  0  Q. 


(5-19) 


where  Q  =  diag  (qfl :  lx  l  weights  the  LOS-error  outputs  that  are  un- 
L  v  u' 

A  CL 

available  for  feedback,  =  diag  (q*) :  NxN  weights  the  displacement 

states  and  Qv  =  diag  (qv] :  NxN  weights  the  velocity  states,  both  assumed 
available  for  feedback  for  the  purpose  of  this  design.  The  control 
weighting  matrix  R:  mxm  in  Eq.  (5-18)  is: 


R  =  (*TBA)Tdiag 


a'' 


(5-20) 


In  general,  R  is  a  fully  populated  matrix. 


Motivation  for  the  use  of  the  cost  functional  (5-18)  incorporating 
an  augmentation  of  the  state  with  LOS-error  outputs  is  connected  with 
the  stringent  tolerance  on  those  outputs.  This  is  discussed  below  in 
Section  5.5.3. 


5.3.2  Application  in  Design 

The  relationships  developed  above  and  in  the  preceding  subsection 

provide  the  foundation  for  a  rational  systematic  approach  to  the 

selection  of  weighting  matrices  for  linear-quadratic  design  in  the 

presence  of  broad-band  stochastic  disturbance  inputs.  The  first  step  is 

to  determine  a  set  of  closed-loop  modal  parameters  wa<j,  £a(j  for  each 

mode  to  be  controlled  based  on  considerations  of  attenuating  the  effects 

of  the  disturbance  input.  Such  determinations  are  discussed  in  Section 

5.2.  Next,  these  desired  values  are  used  with  Eqs.  (5-16)  and  (5-17) 
c  c 

(in  the  place  of  u>  and  C  )  to  determine  the  elements  of  the  matrices 
a  a 

Qd  and  Qy  in  Eq.  (5-19)  as  a  function  of  the  free  parameters  ra.  For 
the  designs  reported  on  here,  parameters  ra  were  chosen  as  (cf.  Eq. 
(5-7)): 
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which  includes  the  average  effect  of  the  disturbance  upon  the  o-th  mode 
and  a  measure  of  the  modal  actuator  influence. 
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Calculation  of  the  parameters  q”,  q”,  and  r  which  generate  the 

D  v  a 

weighting  matrices  are  shown  in  Table  5-2  corresponding  to  the 

closed-loop  modal  parameter  objectives  listed  for  =  5.0 

in  Table  5-1 .  The  LOS-error  weighting  parameters  were  determined 

L 

independently  so  as  not  to  dominate  the  displacement  weighting 

associated  with  the  lowest-frequency  mode  in  the  design  model  (mode  7). 

For  a  single  mode  subsystem  represented  by  Eqs.  (5-13),  any  choice 

of  weighting  matrices  for  the  functional  in  Eq.  (5-14)  satisfying  Eqs. 

(5-16)  and  (5-17)  ensures  that  the  closed-loop  system  formed  by  the 

feedback  law  Eq.  (5-15)  will  have  the  desired  modal  parameters 
c  c 

a>  and  t  used  in  Eqs.  (5-16)  and  (5-17).  However,  since  the  above 
a  a 

analysis  does  not  account  for  the  effect  of  coupling  between  modes  when 
full-state  feedback  is  implemented  on  the  N-mode  design  model,  the 
weighting  matrix  determination  is  not  completely  effective  in  achieving 
the  desired  closed-loop  modal  parameters.  Nevertheless,  it  is  quite 
helpful  in  generating  a  starting  point  for  the  design. 


5.4  Design  Results 

Numerical  experience  substantially  validates  the  selection 
procedure  outlined  In  Section  5.3.  Plots  of  spectral  density  of  the 
LOS-error  outputs  for  the  final  design  run  are  shown  in  Figures  5-1  a 
through  5-1 f.  The  RMS-error  values  associated  with  each  plot,  e.g.. 


RMS  (Wx>  = 


5  (hf 


■J-  /  s  37 

7  IT  *  V 


U)  dot 


»  LOS-X 


are  listed  in  Table  5-3  (last  column).  The  final  closed-loop  pole 
locations,  along  with  associated  damping  ratio  and  characteristic 
frequency,  are  listed  in  Table  5-4.  The  weighting  matrices  Q  and  R  used 
for  the  final  run  are  listed  in  Table  5-5.  (The  column  order  reflects 
the  ordering  of  the  inodes  shown  in  Table  5-1.)  The  final  matrix  Q 
includes  adjustments  in  the  weighting  associated  with  modes  7  and  9  so 
as  to  keep  the  closed-loop  characteristic  frequencies  associated  with 
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these  inodes  above  the  natural  frequency  of  the  isolator  springs.  The 
optimal  state -feedback  gain  matrix  is  listed  in  Table  5-6. 


Numerical  Experience  in  the  Design  Process 


Although  the  properties  of  the  control  system  generated  by  the 
final  design  run  are  of  principal  interest,  some  observations  made  while 
proceeding  toward  the  final  design  are  worth  mentioning.  The  principal 
observations  are  outlined  below. 


5. 5. 1  Optimal  vs.  Nonoptimal  Design 


The  intitial  design  run  incorporates  weighting  on  velocity  states 
only  (i.e.,  Qt  =  0,  Qd  =  0  in  Eq.  (5-19))  using  appropriately  scaled 
values  for  the  parameters  listed  for  wa(3  =  <*iao»  Pad  =  5.0,  *-n 
Table  5-2.  It  turns  out  that  the  control  weighting  matrix  R  generated 
with  these  parameters  by  Eq.  (5-20)  is  ill-conditioned  with  respect  to 
inversion.  This  results  in  failure  to  obtain  a  solution  for  the  Riccati 
equation  associated  with  the  optimization  of  criterion  Eq.  (5-18). 
Nevertheless,  a  stable  design  (relative  to  the  design  model)  is  obtained 
which  gives  surprisingly  close  agreement  with  the  damping  objectives 
calculated  for  the  corresponding  data  in  Table  5-1  (seven  of  ten  pole 
pairs  in  close  agreement).  None  of  the  subsequent  design  runs  produces 
such  close  agreement.  Although  the  RMS  LOS-error  values  for  this  first 
run  are  quite  unsatisfactory,  this  observation  suggests  that  other 
design  techniques  (e.g. ,  direct  pole  assignment)  could  be  worth 
investigating. 


5.5.2  Conditioning  of  the  Control 


Matrix 


Satisfactory  solution  of  the  Riccati  equation  asociated  with  Eq. 
(5-18)  is  achieved  on  subsequent  design  runs  by  using  a  well-conditioned 


control  weighting  matrix  R  obtained  from  the  matrix  of  Eq.  (5-20)  by 
adding  a  scalar  factor  to  the  diagonal  elements: 


r.  .  «-  r .  .  { 1  +a) 
11  11 


This  modification  leads  to  some  deterioration  in  the  agreement  with  the 
modal  parameter  objectives  of  Table  5-1.  The  final  value  of  a  *  0.1 
gives  an  acceptable  level  of  conditioning  to  the  matrix  R  as  well  as 
reasonable  agreement  with  the  modal  parameter  objectives. 


5.5.3  State  vs.  Auqmented-State  Weiqhtin 


No  substantial  reduction  in  the  RMS  LOS -error  values  is  obtained 
until  explicit  weighting  on  those  outputs  is  incorporated  into  the  cost 
functional  in  Eq.  (5-18)  (i.e.,  *  0  in  Eq.  (5-19)).  This  is 

immediately  evident  from  a  comparison  of  the  columns  in  Table  5-3.  Thi 
emphasizes  the  importance  of  formulating  the  design  problem  as  an 
augmented-state  regulator,  in  which  weighting  on  the  specific  linear 
combinations  of  the  state  variables  that  are  of  principal  interest  is 
employed  in  addition  to  the  usual  weighting  on  the  state  variables 
themselves.  This  is  to  be  distinguished  from  an  output  regulator  in 
which  weighting  is  placed  only  on  the  outputs  of  interest,  and  not 
necessarily  on  all  of  the  state  variables  (cf.  Section  4,  Eq.  (4-9)). 


5.5.4  Active  Dampinq  vs.  Active  Stiff ness -Dampinq  Auqmentation 

As  discussed  in  Section  5.2,  the  augmentation  of  damping  by  active 
control  beyond  the  point  of  critical  damping  is  likely  to  be  counter¬ 
productive.  This  constitutes  a  fundamental  constraint  in  the  design  of 
vibration  controllers.  Often,  as  in  the  present  design  problem,  the 
incorporation  of  subcritical  damping  alone  is  not  sufficient  to  meet  the 
design  specifications.  Active  stiffening  in  concert  with  active  d  .mping 
can  overcome  this  constraint  and  lead  to  an  acceptable  design.  In  the 
case  of  a  design  directed  toward  attenuating  the  effects  of  a  broad-band 


disturbance  input,  such  stiffening  has  the  effect  of  moving  the  poles 
associated  with  the  most  disturbance-sensitive  modes  outside  of  the 
disturbance  bandwidth  (cf.  Table  5-1,  last  three  columns).  This  is 
completely  consistent  with  the  corresponding  observation  noted  in  the 
disturbance -rejection  control  design  of  Section  4. 


5.6  Work  in  Progress 

Questions  of  substantial  importance  have  not  been  discussed  in 
this  section.  One  is  the  question  of  stability  of  the  reduced-order 
controller  when  connected  to  evaluation  models  of  larger  order.  Another 
is  the  question  of  the  effect  upon  the  stability  and  performance  of  the 
system  of  incorporating  a  state  observer  into  the  controller.  Both 
questions  are  under  investigation  and  will  be  discussed  in  subsequent 
reports. 
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Table  5-1.  Calculation  of  desired  closed-loop  modal  parameter  objectives 


9.74569 

7.28061 
0. 71608 

3.50193 

3. 74593 

3.86272 

3.99818 
0.93994 

41.2344 

50.6740 


0.001 

0.001 

0.001 

0.001 

0.001 

0.001 

0.001 

0.001 

0.001 

0.001 


1746.220 

513.582 

403.338 

202.140 

108.935 

83.232 

76.559 

50.996 

30.690 

25.829 


9.74569 

7.28061 

0.71608 

3.50193 

3.74593 

3.86272 

3.99818 

0.93994 

41.2344 

50.6740 


1467.56 

488.72 

1184.17 

270.21 

139.99 

104.92 

94.79 

123.81 

19.27 

16.72 


°a0'  Pad 


=  5.0 


C  .  <  0.707;  p 
ad  —  ad 


9.74569 

7.28061 
0.71608 

3.50193 

3.74593 

3.86272 

3.99818 
0.93994 

41 . 2344 

50.6740 


c«d 

pad 

CO  . 
ad 

Cad 

13.8275 

5.0 

35.59 

0.707 

3.89206 

5.0 

14.77 

0.707 

5.25117 

5.0 

1.468 

0.707 

1.27616 

5.0 

4.357 

0.707 

0.43109 

5.0 

3.74593 

0.43109 

0.26071 

5.0 

3.86272 

0.26071 

0. 22215 

5.0 

3.99818 

0.22215 

0. 10399 

5.0 

0.93994 

0.10399 

0.03643 

5.0 

41.2344 

0.03643 

0.02610 

5.0 

50.6740 

0.02610 

Table  5-2.  Quadratic  weights  to  achieve  desired 

closed-loop  modal  parameter  objectives. 


u>  ,  =  < 

ad 

^aO'  pad 

= 

5.0 

0.707;  p 

ad 

= 

5.0 

MODE 

r 

a 

a 

qD 

a 

qv 

a 

qD 

a 

qv 

23 

7.5162 

X 

101 

0 

5.46 

X 

io6 

1  .20 

X 

io8 

1  .43 

X 

104 

22 

8.7760 

X 

io1 

0 

2.819 

X 

io5 

3.93 

X 

106 

9.31 

X 

io3 

7 

1  .9080 

X 

10° 

0 

1 .079 

X 

io2 

8.36 

X 

10° 

1  .96 

X 

10° 

8.4371 

X 

io“2 

0 

6.7403 

X 

10° 

1  .78 

X 

io1 

2.07 

X 

10° 

1 .0980 

X 

io"2 

0 

1.145 

X 

io-1 

0 

1.145 

X 

io-1 

3.8073 

X 

io"4 

0 

1  .5445 

X 

io-3 

0 

1 .5445 

X 

io-3 

15 

8.9494 

X 

io"3 

0 

2.824 

X 

io-2 

0 

2.824 

X 

io-2 

9 

3.1965 

X 

io-5 

0 

1.2215 

X 

io-6 

0 

1  .2215 

X 

io-6 

33 

7.9410 

X 

io'4 

0 

7.162 

X 

io"3 

0 

7.162 

X 

io-3 

34 

3.4933 

X 

io*4 

0 

2.441 

X 

io-3 

0 

2.441 

X 

io-3 

RMS  Error 

( n-raHi  \ 

“ad  ■  “a0f  »ad  *  S’° 

^<  0*7071 

1  »ad  "  S'° 

V0'  Qd=0 

V0'  Qd*o 

v0'  v° 

INPUT 

LOS-X 

87.8635 

89.0878 

0.599137 

AT 

LOS-Y 

282.896 

290.835 

12.7934 

NODE  37 

LOS-Z 

178.351 

165.763 

40.9001 

INPUT 

LOS-X 

139.392 

49.0196 

1 . 33162 

AT 

LOS-Y 

376.895 

174.084 

25.3192 

NODE  46 


LOS-Z 


339.485 


104.179 


79.9046 


Table  5-4.  Closed-loop  poles  for  final  design  run 


NUMBER 

REAL 

IMAG 

DAMPING 

FREQUENCY 

(Rad/Sec ) 

2 

-8. 769678834E-01 

1.461 51 2935E+00 

5.1 4521 691 9E-01 

1.704433257E+00 

4 

-1 .087287502E-01 

3. 408435863E+00 

3. 1 88367784E-02 

3. 41 01 69640E+00 

6 

-1 .948704685E-02 

3. 950036962E+00 

4. 933323384E-03 

3. 950085030E+00 

8 

-3. 8692251 21E+00 

2. 685924262E+00 

8. 21 4731 632E-01 

4. 71 01 05326E+00 

10 

-1 . 089833479E+00 

4. 81 3085291E+00 

2. 208407490E-01 

4.934929283E+00 

12 

-4. 3731 51 345E+00 

5.  428774954E+00 

6. 273270603E-01 

6. 971 086729E+00 

14 

-9. 748363547E+00 

1 . 021 557284E+01 

6. 903695724E-01 

1.41 2050000E+01 

16 

-2. 495696740E+01 

2. 51 0249821 E+01 

7. 050481 302E-01 

3. 53975371 8E+01 

18 

-2. 906650868E+01 

5. 839874058E+01 

4. 455834609E-01 

6. 523246759E+01 

20 

-5. 089422278E+01 

6. 547852743E+01 

6. 1 36887361E-01 

8. 2931 65539E+01 
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N  K>  > 


4.67  x  10®  4.67  x  10® 
1.2  x  10®  3.93  x  106 
1.43  x  104  9.31  x  103 


2.26  x  102  1.78  x  101  0 

1.96  x  10°  2.07  x  10°  1.145  x  10-1 


3.77  x  10 


1.5445  x  10-3  2.824  x  10“2 


7.162  x  10~3  2.441  x  1C)"3 


1 .6204210000-06 
5.293I20000D-06 
1.2603300000-05 
6.1558600000-06 
2.2522900000-06 
2.3473300000-06 
8.7530200000-07 
3.3325500000-06 
3.3828700000-06 


2.3473300000-06 

-1.3470100000-05 

3.4016700000-05 

1.702410000D-05 

5.3564500000-06 

6.2551500000-06 

2.0516900000-06 

3.0029400000-06 

3.0006000000-06 


-5.2931200000- 

3.5751870000- 

-8.2622200000- 

4.143380000D- 

1.2897200000- 

-1.3470100000- 

-4.9140100000- 

5.7894600000- 

-5.7956600000- 


06  1 . 260330000D-05 
05  -8.2622200000-05 
05  2 . 323739000D-04 
05  -1.0614300000-04 
05  -3. 254450000D-05 
05  3. 40167000 0D-05 
06  1.2377700000-05 
06  -1 . 129750000D-05 
06  1.131790000D-05 


8.7530200000-07 

-4.9140100000-06 

1.2377700000-05 

-6.1963100900-06 

-1.9840900000-06 

2.0516900000-06 

8.3222370000-07 

-1.1609500000-06 

1.1716200000-06 


-3.3825500000-06 
5.7894600000-06 
-1 . 129750000D-05 
5.0006400000-06 
2.888830000D-06 
-3.0029400000-06 
-1.1609500000-06 
2.4331010000-05 
-2.0100800000-05 


-6.1558600000-06 

4.1433800000-05 

-1.0614300000-04 

5.8739120000-05 

1.630920000D-05 

-1.7024100000-05 

-6.1963100000-06 

5.0006400000-06 

-5.0166400000-06 


3.3828700000-06 

-5.7956600000-06 

1.1317900000-05 

-5.0166490000-06 

-2.8948300000-06 

3.0006000000-06 

1.1716200000-06 

-2.0100800000-05 

2.4360270000-05 


-2.2522900000-06 
1.2897200000-05 
-3.2544500000-05 
1.6309200000-05 
5 . 734795000D-06 
-5.3564500000-06 
-1.9840900000-06 
2.8838300000-06 
-2.8948300000-06 
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Table  5-6.  Optimal  state-feedback  gain  for  final  design  run 


1 

1 . 1443134640*06 
-3. 4149455910*05 
1.4627233130*05 
-2.4641875850*05 
-3.8651743100*05 
5.1599628100*05 
1.1294570800*06 
-1.7819657110*04 
1.4181279770*04 

6 

1.2860790660*04 
-1. 5856836 35D+03 
8.6125703960*02 
-5.9036577590*02 
4.7379915800*03 
-6.5777229000*03 
-1.1965549740*04 
9.6934034220*02 
-9.0993254450*02 

11 

5.7232685720*04 

-1.6477178850*04 

7.2806852190*03 

-1.1568380910*04 

-1.1629267380*04 

1.601101745D*04 

3.5763766060*04 

-5.0760603090*02 

4.4721400480*02 

16 

2.6507608100*02 

-3.4401666360*01 

1.1461465230*01 

-1.4522821540*01 

7.3482834620*01 

-1.1656100610*02 

-2.9874664870*02 

5.8046020060*00 

-3.0097798390*01 


2 

1.4073392790*06 

3.7176301830*04 

•2.6882831150*04 

7.3673401940*04 

6.1377414390*04 

-7.1881842490*04 

-6.8334786110*04 

-3.5833418390*04 

3.2136976460*04 

7 

-5.0364405430*02 

1.1122719740*02 

-6.1207844790*01 

2.8487732730*01 

-3.5682433870*03 

-4.4043424030*03 

5.3216697930*03 

1.1595644130*02 

1.1079828870*02 

12 

1.9809579320*05 
6.3070190060*03 
-4.4427842770*03 
1 . 1012345330*04 
5.1985235140*03 
-5.4821611190*03 
-5.6850547330*02 
-2.5281004850*03 
2.4361327200*03 

17 

5.8847632400*00 

2.0225441950*01 

1.0208794430*00 

-2.4218692210*00 

-3.8383152740*01 

-5.4753631600*01 

1.4480026030*02 

-1.0538686330*02 

-1.0917901280*02 


3 

-1.4167218410*04 

5.0704091400*02 

-2.0066634670*02 

-1.8252114240*02 

•1.9128564160*03 

2.5677529770*03 

4.648471786D*03 

-1.7133243880*03 

1.6501576680*03 

8 

-3.2848560470*02 

7.5734189470*01 

-4.1689055150*01 

1.6946331670*01 

-2.5441203990*03 

-3.2100458280*03 

3.6894746570*03 

-2.7227918000*02 

-2.8175135910*02 

13 

-1.1919342630*04 

-1.8362699610*02 

1.8191165190*02 

-5.4537697270*02 

-3.3366473730*01 

2.7124672970*01 

-3.5885383030*02 

-1.8320116430*03 

1.5776214600*03 

18 

-8.628057114D*00 

1.1539021430*00 

-1.085416906D*00 

1.025622170D*00 

-2.5253188530*01 

-5.0897271340*01 

1.1844868700*02 

-1.306291251D*02 

-1.3108695310*02 


4 

-1 . 733328386D*04 
2.1975713650*03 
-1.2007565970*03 
6.5444171700*02 
-6.654043762D*03 
9.0012776700*03 
1.6807154550*04 
-1.3939409600*03 
1.0216693000*03 

9 

3.3788619380*04 

-7.6752930080*03 

4.0729829210*03 

-1.5726728930*03 

2.6971542580*05 

3.3556271470*05 

-3.8659425380*05 

1.1734156060*03 

1.6759470360*03 

14 

1.6663982210*03 
1 . 5044704200+02 
-4.8708133290*01 
1.3807316680*02 
-5.403490472D*01 
2.1579657520*02 
1.1658815550*03 
2.9478034270*02 
-2.5062017410*02 

19 

6.3857568660*02 

-1.4173648610*02 

7.4654604750*01 

-2.606454526D*01 

5.2715750710*03 

6.6328705220*03 

-7.5196665350*03 

2.1457579940*01 

3.0029622510*01 


5 

-1.6522672150*03 
3.9083912120*02 
-2.1400647260*02 
9. 1317452360*01 
-1.3560314670*04 
-1.6846148930*04 
1.9589210230*04 
8.042224918D*02 
7.5373747820*02 

10 

7.1848824870*04 

-2.4564261660*04 

1.4019376150*04 

-1.3220540210*04 

5.609774694D*04 

-1.0801006900*05 

-1.6688068930*05 

-1.4262920430*03 

1.7987830940*03 

15 

-5 . 445847648D+01 
-9.6911345870*00 
-8.1339241550*00 
9.1709432420*00 
-1.3266486760*02 
-2.846416115D+02 
6.459121099D*02 
-3. 756881147D+02 
-3.7665581640*02 

20 

2.777572758D+03 

-6.5786362380*02 

4.8621050750*02 

-4.5529481670*02 

2.0640461500*03 

-3.6098367830*03 

-5.7976812500*03 

-2.9394325570*01 

3.5429014400*01 
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Figures  5-1.  Output  power  spectral  density  from  closed- 
loop  linear  quadratic  regulator 
(augmented-state  formulation). 
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igure  5-la.  Disturbance  at  node  37,  output  LOS-: 


SECTION  6 


LINEAR-QUADRATIC  OPTIMAL  CONTROL  DESIGN 
WITH  SPILLOVER  REDUCTION 


6.1  Introduction 

Both  Sections  4  and  5  describe  the  design  of  optimal  control 
systems  for  the  10-mode  model  composed  of  the  top  10  with  respect  to 
their  RMS  LOS-errors  (Modes  23,  22,  7,  12,  13,  14,  15,  9,  33,  and  34  in 
the  descending  order).  Spillover  (and  possibly  closed-loop  instability) 
is  expected  when  these  control  systems  are  applied  to  the  20-mode 
evaluation  (or  "truth")  model,  which  contains,  in  addition,  the  next  top 
10  modes  (modes  45,  40,  37,  16,  24,  10,  47,  48,  62,  and  56,  again  ir.  the 
descending  order).  In  this  section,  we  report  a  design  of  actuator 
synthesizers  for  preventing  spillover  of  selected  low-frequency  modes 
and  a  corresponding  linear-quadratic  optimal  control  design. 


A  Design  of  Actuator  Synthesizers 


To  demonstrate  the  design  of  a  spillover-preventing  synthesizer 
for  the  given  six  actuators,  we  consider  modes  23,  22,  7,  12,  13,  14, 

15,  9,  33,  and  34  as  the  "primary  modes"  and  treat  modes  10,  16,  24,  37, 
40,  and  45  as  the  "secondary  modes".  In  fact,  all  are  critical  to  the 
LOS  performance  and  should  be  actively  controlled  if  the  design  model  is 
not  limited  to  only  10  modes.  Note  that  the  above  secondary  modes  are 
of  the  lowest  frequencies  among  the  10  nonprimary  inodes  and  that  inodes 
10,  16,  and  24  lie  between  the  closely-spaced  primary  modes.  The 
actuator  synthesizer  is  designed  to  achieve  two  objectives:  (1)  to 


I 

i 
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prevent  control  spillover  to  the  six  secondary  modes,  and  (2)  to  correct 
the  modal  control  authority  of  the  four  most  LOS-critical  modes,  namely, 
modes  23,  22,  7,  and  12.  Both  objectives  are  achieved  :  tisfactorily  by 
a  synthesizer  r  which  is  decomposed  into  two  parts  as 

T  =  Q2T  (6-1) 

The  matrix  Q2  is  computed  first  (for  the  first  objective)  followed  by 

* 

the  computation  of  r  (for  the  second  objective). 

6.2.1  Synthesis  for  Prevention  of  Control  Spillover 

Ideally,  we  wish  to  synthesize  to  prevent  control  spillover  to  the 

six  secondary  modes  exactly,  but  we  decided  to  compromise  with  the 

control  of  the  top  10  LOS-critical  modes.  Since  the  6x9  actuator 
T 

influence  matrix  $gBF  corresponding  to  the  six  secondary  modes  has  rank 

6,  an  exact  prevention  of  control  spillover  to  the  six  secondary  modes 
would  require  6  independent  actuators  and  hence  result  in  at  most  3 
independent  synthetic  actuators  (see  Ref.  6-1).  Note  that  all  the  nine 
actuators  are  used  to  control  the  10  primary  modes  when  no  prevention  of 
spillover  is  considered.  We  felt  that  it  would  be  more  desirable  if  one 
more  independent  synthetic  actuator  was  reserved  for  controlling  the 
primary  modes.  With  such  a  restriction,  the  synthesizer  could  be 
designed  in  one  of  the  following  alternative  ways;  (1)  ignore  one 
secondary  mode  and  prevent  control  spillover  to  the  remaining  5 
secondary  modes  exactly,  or  (2)  prevent  spillover  to  all  the  six 
secondary  modes  approximately.  We  took  the  second  alternative,  since  a 
"degree-of-independence"  test  revealed  that  the  actuator  influences  on 
these  secondary  modes  have  a  high  degree  of  dependence.  So,  some 
spillover  to  each  of  the  six  secondary  modes  should  be  expected. 


The  exact  synthesis  program  previously  reported  in  Ref.  6-2  was 
modified,  and  a  9  x  4  matrix  Q2  that  provides  approximate  spillover 


prevention 

thus  computed. 

The  following 

shows  the  result: 

0.40377D-01 

0.68237D+00 

0.0 

0.0 

-0.21 362D+00 

0.4201 7D+00 

-0 . 70978D+00 

0. 35879D+00 

0.71674D+00 

0.78608D-01 

0.11 774D+00 

0.45053D+00 

0.1  3984D+00 

0.82051D-01 

0.81981D-01 

0. 39338D-01 

Q2  - 

-0.91406D-01 

0.38219D+00 

0 . 29262D+00 

-0.31 405D+00 

0. 54240D-01 

-0.40540D+00 

-0 . 3001 8D+00 

0.29277D+00 

-0 . 37592D+00 

-0.1 8267D+00 

-0.17103D+00 

-0.1 4506D+00 

-0.3651 3D+00 

0.23773D-01 

0.36782D+00 

0.48031 D+00 

_ -0 . 36542D+00 

0 .237370-01 

0.36790D+00 

0.48028D+00  _ 

6.2.2  Synthesis  for  Correction  of  Modal  Control  Authority 

Modes  23,  22,  7,  and  12  are  the  most  LOS-critical  primary  modes: 
their  individual  contributions  to  LOS  errors  not  only  rank  the  highest 
among  the  20  modes  (also  among  all  the  156  modes  of  the  VCOSS  model 
structure),  but  are  also  disproportionately  large  compared  with  other 
LOS-critical  modes.  For  example,  they  are,  respectively,  34,  10,  8,  and 
4  times  larger  than  that  of  Mode  9.  Therefore,  the  control  authority  of 
the  9  actuators  on  these  four  modes  not  only  should  rank  the  highest  but 
also  should  have  the  similar  proportions  as  compared  with  that  on  mode 
9,  for  example.  By  whatever  measure  of  control  authority,  the 
situations  are  just  the  opposite,  however.  As  can  be  seen  from  Tables 
6-1  through  6-4,  these  four  modes  all  have  very  disappointingly  small 
values,  and  all  lie  at  or  near  the  bottom,  even  with  the  nonprimary 
modes  included  for  comparison. 

A 

The  free  parameter  matrix  T  is  useful  for  trying  to  correct  the 

A 

situation.  We  calculated  the  matrix  r  column  by  column  using  a  least- 
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p, . 
’  V 


squares  optimization  technique  in  such  a  way  that  each  column  tried  to 


i 

T7- 

i-.-. 


elevate  the 

relative  control  authority 

cm  one  of  these  four  modes. 

The 

results  are 

shown  below: 

'  ■ 

-0.58244D+01 

0.971 29D+00 

-0.24545D+00 

0.11 582D+00 

• 

A 

0. 23931 D+02 

-0.451 18D+01 

0.53743D+00 

-0.90634D+00 

p 

r  = 

-0.19278D+02 

0.34960D+01 

-0.52266D+00 

0.60541 D+00 

}•. 

0.91 904D+01 

-0.1 6749D+01 

0. 22756D+00 

-0 . 29022D+00  _ 

i 

Forming  the  matrix 

product  of  Q2 

and  f  as  in  Eg 

.  (6-1  ) ,  we 

obtained  the  synthesizer 

T: 

0.1 6094D+02 

-0.30395D+01 

0.35681 D+00 

-0.61 378D+00 

9 

0.28280D+02 

-0.51 856D+01 

0.73087D+00 

-0.93940D+00 

-0.42270D+00 

-0.1 4559D-02 

-0.92692D-01 

-0.47700D-01 

i- 

-0.69848D-01 

-0.1 365 ID-01 

-0.241 23D-01 

-0.1 9953D-01 

\!  ’ 

r  = 

0.11511 D+01 

-0.2641 5D+00 

0.34306D-02 

-0.88683D-01 

1 

-0.15400D+01 

0. 34201 D+00 

-0.76766D-02 

0.1 0702D+00 

-0.21 796D+00 

0.1 0408D+00 

0.50477D-01 

0.60575D-01 

0.1 9005D-01 

0.19527D-01 

0.1 9451 D-01 

0.19446D-01 

•  ■ 

_  0.1 7885D-01 

0.19762D-01 

0.19451D-01 

0.19506D-01 _ 

■ 

,V 

The  synthesized  actuator  influences  are  thus  given  by  post- 

multiplying 

the  unsynthesized  influence 

matrix  by  the 

matrix  T.  Tables 

6-5  through 

6-8  show  the 

ranking  of  the 

synthesized  control  authority 

’  • 

* 

using  the  same  measures  as  before.  Observe  that  the 

synthesized  control 

i  ‘ 

authority  on  the  two  most  LOS-critical 

modes,  namely 

23  and  22,  have 

been  greatly  improved,  compared  with  that  on  mode  9,  for  instance, 
though  improvement  in  modes  12  and  7  are  rather  small.  The  proportions 
of  the  synthesized  control  authority  among  the  primary  modes  are  also 
improved . 


n 


The  difficulties  in  improving  the  control  authority  on  modes  12 
and  7  manifest  two  levels  of  improper  placement  of  the  actuators*  The 
first  level,  shown  even  clearer  by  modes  23  and  22,  is  the  undesirably 
small  magnitude  of  actuator  influences  on  the  most  important  modes.  The 
second  level  is  the  undesirably  high  dependence  of  the  primary-mode 
actuator  influences  on  the  secondary- mode  influences. 

6.3  Optimal  Control  Design  with  Spillover  Reduction 

With  the  spillover-preventing  synthesizer  r  included,  the  dynamic 
equations  become 

\  +  2?i“i\  +  “ini  =  ♦iBFru  (6“2) 


for  each  primary  mode.  The  only  difference  is  the  replacement  of  the 
unsynthesized  influence  matrix  Bp  by  the  synthesized  matrix  Bpr. 

Thus,  any  controller  design  technique  that  is  applicable  to  the 
unsynthesized  case  is  equally  applicable  to  the  synthesized  case.  In 
particular,  both  the  disturbance-rejection  design  presented  in  Section  4 
and  the  linear-quadratic  optimal  control  design  presented  in  Section  5 
can  be  mimicked  here  for  the  synthesized  case.  We  chose  to  mimic  the 
latter  in  our  first  attempt  of  reduced-order  modern-controller  design 
with  spillover -preventing  synthesizers.  The  analytical  and  one-shot 
nature  is  both  cost-effective  and  time-effective,  and  hence  is  more 
attractive  than  the  regular  trial-and-error  approach.  It  is  based  on 
the  analytical  procedure,  jointly  developed  with  J.G.  Herther  earlier 
[Ref.  6-3],  for  selecting  the  weight  matrices  in  the  linear-quadratic 
regulator  design. 


ladratic 


Like  the  design  of  Section  5,  the  objective  is  to  attenuate  the 
individual  contribution  to  the  RMS  LOS-error  from  each  primary  node  to  5 
yrad .  Note  that  the  defocus  error  (i.e.f  the  "LOS-Z"  error)  has  not 
been  used  in  the  ranking  and  selection  of  the  10  control-system  design 
model  modes,  i.e.,  the  10  primary  modes.  Note  also  that  the  total 
contribution  to  the  RMS  LOS-error  from  all  the  10  primary  modes  is 


/l 0  x  52  =  15.81  yrad 

if  all  the  individual  contributions  are  reduced  to  5  yrad. 

Following  the  same  analytical  procedure  as  sketched  in  Section  5 
(see  Eqs.  (5-18)  to  (5-21),  and  refer  to  Eqs.  (5-11),  (5-12),  (5-16), 
(5-17),  and  the  case  of  £a<j  £  0.707  in  Table  5-1),  we  selected  a 
priori  the  weight  matrices  Q  and  R  as  shown  in  Table  6-9. 

6.3.2  Design  Results 

Table  6-10  shows  the  resulting  state-feedback  gains.  Table  6-11 
lists  the  resulting  closed-loop  poles,  modal  damping  and  modal 
frequencies.  The  power  spectral  density  of  the  LOS-errors  subject  to 
the  broad-band  random  disturbances  at  nodes  37  and  46  are  shown  in 
Figures  6-1  a  through  6-1 f.  Table  6-12  summarizes  the  RMS  LOS-errors. 

Observe  that  the  total  LOS-Y  error  contributed  from  all  the  10 
modes  subject  to  either  or  both  disturbances  is  very  close  to  the  design 
objective  of  5  yrad  for  each  mode.  Although  the  total  LOS-X  error  has 
also  been  reduced,  it  requires  further  reduction  to  meet  the  design 
objective.  Some  analysis  and  adjustment  of  the  LOS  weighting  are 
required  after  this  preliminary  design. 

Recall  that  in  this  design,  only  four  control  inputs  are  used  for 
controlling  the  10  primary  modes,  in  contrast  with  the  designs  presented 
in  Sections  4  and  5,  in  which  all  nine  control  inputs  are  used.  Lower 
performance  should  be  expected. 


•nil 


Other  detailed  discussions  as  well  as  the  follow-on  design 
activities  will  be  reported  later. 
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Table  6-3.  The  damped  degree-of-controllability  (DOC)  values  for 
each  mode's  control  influence  with  delta  *  .06. 


RANK 

MODE 

DOC 

1 

9 

0.84537191 88D-05 

2 

10 

0.5165528493D-05 

3 

14 

0.1 8846 3244 3D-05 

4 

16 

0 . 4559450893D-06 

5 

45 

0 . 2046270947D-06 

6 

48 

0.1 975475353D-06 

7 

34 

0.11 33468352D-06 

8 

13 

0.111 9594247D-06 

9 

56 

0.111 6967443D-06 

10 

33 

0.8376070304D-07 

1 1 

47 

0.69331 73299D-07 

12 

15 

0 . 6781 903938D-07 

13 

40 

0.6642650620D-07 

14 

62 

0.544911 21 99D-07 

15 

12 

0.501 9237498D-07 

16 

37 

0. 35971 69420D-07 

17 

7 

0 . 8860244364D-08 

18 

24 

0. 5068259371 D-08 

19 

23 

0.41 23635709D-08 

20 

22 

0. 30853561 73D-09 

Table  6-4.  The  undamped  DOC  values  for  each  mode's  control 
influence  with  delta  *  1 . 


RANK 

MODE 

DOC 

1 

9 

0.221 4876900D-02 

2 

10 

0.1 324474663D-02 

3 

14 

0 . 2320624786D-03 

4 

16 

0 . 5358790490D-04 

5 

13 

0.1 42583301 5D-04 

6 

15 

0.8042876959D-05 

7 

12 

0.6879592027D-05 

8 

45 

0.347801521 6D-05 

9 

48 

0 . 3291 040802D-05 

10 

7 

0.23788651 93D-05 

11 

34 

0.1 88435781 2D-05 

12 

56 

0.1 860633769D-05 

13 

33 

0 . 1 438250407D-05 

14 

47 

0.1 1 62451 603D-05 

15 

40 

0.11 28703537D-05 

16 

62 

0.8998399631  D-06 

17 

37 

0.5991 474092D-06 

18 

24 

0.2244366462D-06 

19 

23 

0. 2080305 574D-06 

20 

22 

0.2050291 577D-07 
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Table  6-5.  The 

RSS  values 

for  each  mode's  control  influence. 

n 

RANK 

MODE 

RSS 

1 

62 

0.7665261 41 2D-01 

2 

34 

0. 1964353669D-01 

3 

56 

0. 1 949282657D-01 

4 

48 

0.1 868951 163D-01 

m 

5 

23 

0. 1 835832973D-01 

P 

6 

47 

0. 588549291 4D-02 

‘ 

7 

9 

0.491 541 361 0D-02 

i 

i 

8 

15 

0.3291 04761 8D-02 

• 

9 

22 

0.291 61 99034D-02 

•'  ■ 

10 

33 

0. 2021 945935D-02 

A 

11 

14 

0. 1  31 001 1 534D-02 

12 

13 

0.11 7581 1702D-02 

13 

12 

0.  8784900923D-03 

' 

-■*  > 

14 

16 

0.  4230897846D-03 

15 

7 

0.231 481 5954D-03 

.  ■' 

16 

10 

0. 1 249775007D-03 

i 

17 

24 

0.  5528453623D-04 

18 

40 

0. 1373992454D-05 

19 

37 

0.  747388731 8D-06 

» 

• 

20 

45 

0. 6335282574D-06 

: 

• 

A 

«  „• 

Table  6-6.  The  terms 

with  the 

greatest  absolute  value  in  each  row. 

i 

-r 

RANK 

MODE 

RSS 

v 

• 

- 

1 

62 

-0. 752028581 9D-01 

■ 

2 

56 

0. 1 91 391 9301D-01 

3 

34 

-0. 1 883667488D-01 

4 

48 

0.1 850321 486D-01 

i, , 

r  . 

5 

23 

0. 180551 6365D-01 

» 

6 

47 

-0. 5652007766D-02 

7 

15 

-0. 3257844598D-02 

»*■  — ■ 

># 

8 

22 

-0. 2857453201D-02 

9 

9 

0. 250491 1298D-02 

m 

s.’  - 

10 

33 

0. 1 996049546D-02 

,1 

/*, 

11 

14 

-0.  1287734647D-02 

,1 

•  ,t 

12 

13 

0. 96805971 92D-03 

» 

13 

12 

0.8641 99361 9D-03 

-i 

14 

16 

-0. 41 68865235D-03 

r* 

15 

7 

-0.  2276252051 D-03 

* 

p.. 

16 

10 

-0. 1231 451 046D-03 

17 

24 

-0.5447396499D-04 

pv 

18 

40 

0. 1353847241D-05 

Lv#- 

19 

37 

-0. 7364306620D-06 

* 

20 

45 

0. 6242395879D-06 

i: 

> 

6-10 

c 

1 

Li- 

’/v  *,*  •- 

*'*I*"^'  in*  '>  — L— J  -a 

Table  6-7.  The  damped  DOC  values 
with  delta  =  .06  . 

RANK  MODE 


1 

62 

2 

23 

3 

34 

4 

48 

5 

56 

6 

9 

7 

47 

8 

15 

9 

22 

10 

33 

1 1 

14 

12 

13 

13 

12 

14 

16 

15 

7 

16 

10 

17 

24 

18 

40 

19 

37 

20 

45 

for  each  mode's  control  influence 


DOC 

0.1 7685091 74D-06 
0.25641 13362D-07 
0. 1 999705975D-07 
0.1 336721 607D-07 
0. 1 322839465D-07 
0.2965294086D-08 
0. 143081 91 56D-08 
0. 1 060278071D-08 
0. 65422201 80D-09 
0. 21 21 223747D-09 
0. 1 839420675D-09 
0.1 51 97261 59D-09 
0. 8737861 122D-10 
0. 1 71 921 681  ID-1  0 
0. 7320702038D-1 1 
0. 1 883760480D-1 1 
0. 2228680891D-1 2 
0. 8565986263D-1 6 
0. 2704895235D-1 6 
0. 1 695593384D-1 6 


Table  6-8.  The  undamped  DOC  values 
with  delta  =  1. 


RANK  MODE 


1 

62 

2 

23 

3 

9 

4 

34 

5 

48 

6 

56 

7 

15 

8 

22 

9 

14 

10 

47 

11 

13 

12 

12 

13 

33 

14 

16 

15 

7 

16 

10 

17 

24 

18 

40 

19 

37 

20 

45 

for  each  mode's  control  influence 


DOC 

0. 1896007763D-04 
0. 1 728207773D-04 
0. 6577580881D-05 
0. 3807248873D-05 
0. 1 84941 0073D-05 
0. 1 41 8432560D-05 
0.1 3300821 78D-05 
0.5801 31 4068D-06 
0. 21 88709486D-06 
0. 1 941 872078D-06 
0. 1 824025239D-06 
0. 1 096285936D-06 
0. 49571 29057D-07 
0. 21 78425043D-07 
0. 1 494975784D-07 
0. 41 61 087073D-08 
0.13662681 80D-09 
0. 1 300632226D-1 3 
0. 493379441 6D-1 4 
0.  26971 99398D-1 4 
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6-9.  Weight  matrices  Q  and  R 
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Table  6-11.  Closed-loop  poles,  modal  damping,  and  modal  frequencies 
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Disturbance 
at  Node  37 
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at  Node  46 
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408.44 

Total 


610.57 


8.48 


610.63 


so.o 


FREQUENCY  (HERTZ) 


0*001- 


6-18 


FREQUENCY  (HERTZ) 


FREQUENCY  IHERTZ1 


0*001 


Figure  6-lf.  Disturbance  at  node  46,  output  LOS- 


SECTION  7 


LARGE-ANGLE  SPACECRAFT  SLEWING  MANEUVERS 

7.1  Introduction 

This  section  presents  some  recent  extensions  for  the  problem  of 
maneuvering  a  flexible  spacecraft  through  large  angles.  In  particular, 
this  section  addresses  specific  control  issues  in  the  following  four 
areas: 

(1)  Selection  of  a  reduced  order  model; 

(2)  Retargeting  maneuvers; 

(3)  Optimal  control  employing  state  and  control  feedback;  and 

(4)  Fixed  time  maneuvers  using  feedback  control. 

The  remainder  of  this  section  is  presented  in  four  parts.  Section 

7.2  presents  a  method  for  selecting  a  reduced  order  plant  model  using  a 
degree-of-controllability  measure  and  a  residual  energy  method.  Section 

7.3  presents  a  technique  for  slewing  a  flexible  spacecraft  to  engage  a 
moving  target.  This  material  represents  an  extension  to  the  work  first 
presented  in  Ref.  7-1.  Section  7.4  presents  a  control  design  using 
integral  precompensation,  where  the  state,  controls,  and  control  rates 
are  fed  through  optimal  gains  into  a  cascade  of  integrators.  The  output 
of  the  last  integrator  is  the  control  which  maneuvers  the  vehicle.  This 
material  represents  an  extension  to  the  feedback  case  of  the  open-loop 
control-rate  penalty  technique  presented  in  Ref.  7-1  . 


Given  the  typically  high  dimension  of  finite  element  models  used 
to  characterize  flexible  spacecraft,  it  is  necessary  to  reduce  the 
flexible-body  plant  description  to  a  manageable  number  of  degrees  of 
freedom.  The  reduced-order  model  selection  procedure  presented  in  this 
section  consists  of:  (1)  computing  a  relative  degree-of-controllability 
for  each  flexible  mode;  and  (2)  using  a  residual  energy  technique  to 
minimize  the  number  of  degrees  of  freedom  retained  in  the  plant 
description. 

For  a  structure  with  many  vibration  modes,  the  relative  degree-of- 
controllability  measure  is  obtained  by  checking  the  rank  of  the  control¬ 
lability  matrix 


(4x4)  (4  x  N  ) 

c 


and  uii  is  the  frequency  of  the  i-th  mode,  V  (2xNc)  is  the  matrix  of 
control  influence  coefficients  for  the  rigid  body  and  the  i-th  mode,  N 
is  the  number  of  modes  in  the  plant  model  before  model  reduction,  and 
Nc  is  the  number  of  independent  controls. 


To  determine  the  rank  of  the  controllability  matrix,  a  Gram- 
Schmidt  orthogonalization  process  is  carried  out  on  the  columns  of  the 
C^.  As  a  result,  the  rank  of  is  found  by  computing  the  number  of 
non-zero  columns.  Numerically,  the  number  of  non-zero  columns  is 
established  by  computing  the  sum  of  the  squares  of  the  coefficients  of 
each  column  and  comparing  the  result  with  a  prescribed  tolerance  e.  By 
varying  the  tolerance  e  and  observing  the  rank  of  as  a  function  of 
e,  one  obtains  a  rough  indication  of  the  relative  degree-of- 
controllability  of  the  i-th  mode.  In  practice  all  modes  which  have  been 
found  to  be  uncontrollable  are  removed  from  the  candidate  set  of  inodes 
considered  for  the  reduced-order  model. 

Using  this  approach  on  an  84-mode  subset  of  the  flCOSS  Model  No.  2 
structure,  20  modes  were  eliminated. 

To  further  reduce  the  number  of  degrees  of  freedom  in  the  plant 
model,  a  residual  energy  method  is  next  presented. 

The  first  step  in  the  residual  energy  method  consists  of 
arbitrarily  selecting  several  of  the  significant  modes  identified  in  the 
degree-of -controllability  measure  as  a  candidate  reduced-order  model. 
Next,  basing  a  control  design  on  the  candidate  reduced-order  model,  the 
resulting  control  is  applied  to  all  the  modes  remaining  after  the 
degree-of -controllability  measure.  At  the  end  of  the  maneuver  the 
control  is  turned  off  and  the  residual  energy  of  the  i-th  mode 

Ei  =  ^  n?(tf)  +  2  “ini(tf)  (i=1»****n  >  (7-2) 

• 

is  computed,  where  is  the  i-th  modal  amplitude,  is  the  time 
derivative  of  the  i-th  modal  amplitude,  is  the  i-th  modal 
frequency,  and  n*  is  the  number  of  modes  remaining  after  the  degree-of- 
controllability  measure.  If  any  uncontrolled  modes  are  observed  to  have 
unacceptably  large  residual  energies,  then  the  candidate  reduced-order 
model  is  expanded  to  include  these  modes. 


As  a  check  on  the  new  reduced-order  model,  a  new  control  design  is 
established,  and  the  resulting  control  is  applied  to  the  structure.  As 
before,  the  energy  in  the  i-th  residual  mode  is  computed  using  Eq. 

(7-2),  the  residual  mode  energies  for  i=1 ,  •••,  n*  are  ranked 
according  to  magnitude,  and  the  reduced-order  plant  model  is  refined  as 
necessary.  Typically,  only  one  or  two  iterations  of  the  residual  energy 
method  are  required  before  a  suitable  reduced-order  model  has  been 
found. 


7.3  Retargeting  Maneuvers  For  Flexible  Spacecraft 

The  problem  of  slewing  a  flexible  spacecraft  to  engage  a  moving 
target  is  presented  in  this  section.  This  material  represents  an 
extension  to  the  rigid  body  retargeting  maneuvers  first  presented  in 
Ref.  7-1  .  In  particular,  the  results  of  this  section  are  valid  for  the 
special  case  of  a  single-axis  maneuver. 

The  target  motion  complicates  the  slewing  control  problem,  because 
the  terminal  control  state  must  be  set  in  order  to  have  the  spacecraft 
track  the  target  motion  at  the  end  of  the  maneuver.  To  obtain  the 
terminal  control,  the  spacecraft's  attitude  motion  is  constrained  to 
osculate  the  target  attitude  motion  through  (k+1 )  time  derivatives, 
where  k  denotes  the  highest  time  derivative  of  the  control  appearing  in 
the  performance  index.  In  addition  to  the  rigid  body  attitude  motion 
constraints,  it  is  convenient  to  impose  the  vibration  constraint  that 
the  flexible  body  response  is  zero  through  (k+1)  time  derivatives. 
Subject  to  the  desired  target  motion  constraints  above,  the  required 
terminal  control  state  is  shown  in  this  section  to  be  a  function  of  the 
target  attitude  motion. 

The  assumed  target  motion  kinematics  and  geometry  is  shown  in 
Figure  7-1 . 


itimal  Control  Problem 


7.3.1  0£ 

The  optimal  control  problem  for  the  retargeting  maneuver  is  to 
seek  the  torque  time  history  ji(t)  which  minimizes  the  performance  index 


1  ,  f  r  T  T 

J  =  -  J.  [s  W  s  +  u  W  u  +  u  W  u  dt 

2  J  t  L—  ss—  —  uu—  —  w— J 


(7-3) 


subject  to 


s  =  As  +  Bu,  given  s(t  )  and  s(t_)  (7-4) 

—  —  —  —  o  —  r 

r  000  -.T  • 

where  =  [0  r  •••nn  0  n1  •  •  *r>n J  »  0  is  the  rigid  body  angle,  0  is  the 
angular  rate,  for  i=1 ,  ••*,  n  denotes  the  generalized  coordinates  for 

the  flexible  degrees  of  freedom,  for  i=1 ,  •••,n  denotes  the  time 

derivative  of  the  generalized  coordinates  for  the  flexible  body  degrees 
of  freedom,  n  is  the  number  of  flexible  body  generalized  coordinates,  ii 
denotes  a  distributed  set  of  independent  controls,  A  is  the  plant 
dynamics  matrix,  B  is  the  control  influence  matrix,  W8S  is  the  state 
weight  matrix,  Wuu  is  the  control  weight  matrix,  and  Ww  is  the 
control-rate  weight  matrix. 

As  shown  in  Ref.  7-1,  the  necessary  conditions  for  the  optimal 
control  problem  follow  as: 


State: 

£ 

=  A  s  +  B  u 

(7-5) 

• 

T 

Costate: 

X 

=  -W  s  -  A  X 

(7-6) 

ss—  — 

Control: 

9 

U 

=  CU  +  DX 

(7-7) 

where 
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Lw_1w 

w  uu 


(2N  x2N  ) 
c  c 
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,  W_1B  , 
L  w  J 


(N  x2N) 
c 


U 


•TiT 

«  J 


Nc  is  the  number  of  independent  controls,  and  N  =  1  +  n. 


The  necessary  conditions  of  Eq.  (7-5),  (7-6),  and  (7-7)  define  a 
two-point  boundary- value  problem;  since  s^  is  known  initially  and 
finally,  _X  is  unknown  initially  and  finally,  and  u^  is  unknown  initially 
and  finally.  However,  as  shown  in  kef.  7-1,  the  terminal  boundary 
conditions  for  ii  can  be  specified  to  suit  the  needs  of  individual 
problems.  In  particular,  in  this  section  the  final  value  of  the  control 
is  set  in  order  to  enforce  the  desired  target  attitude  motion 
constraints  and  the  vibration  suppression  constraints. 

Assuming  that  vi(t0)  is  given,  the  equation  defining  the  solution 
for  ^(t  )  and  _u(t  )  is  given  by  [Ref.  7-1] 


^s  X  ^su. 


Tu  X  u  u, 
a  a  b 


x(to)  ;  \  -ctf) 


u.  (t  ) 
— b  o 


u  (t J' 
—a  f 


$  4* 

rss  su 


4  <j> 

Tu  s  Tu  u 


1  / 

■ 

s(t  )  ) 

r  ° 

1  u  (t  )\ 

J  1 

—a  o  j 

(7-8) 
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,  A  A  • 

where  u  =  u  ,  u  =  u  ,  and  ($.,<(>  ,  d>  .,  ,  d>  ,  <f>  ,  <*>  , 

— a  —  — b  —  rsX  rsu.  yu  X  *u  a  *ss'  *su  '  s' 

b  a  a  b  a  a 

$u  )  denote  blocks  of  the  state  transition  matrix  for  Eqs.  (7-3), 
a  a 

(7-4),  and  (7-5). 

Before  Eq.  (7-8)  can  be  inverted  for  _X(t0)  and  Ub(tQ),  the 
terminal  control,  Ua<tf)»  must  be  specified.  The  terminal  control 
for  the  retargeting  maneuver  is  obtained  from  the  vehicle's  equation  of 
motion. 

In  order  to  illustrate  the  technique,  the  vehicle's  equation  of 
motion  is  assumed  to  be  given  by 


— 0n  9(t) 


— 0n  %  h(t) 


0  K 


1  4M  Vfc) 


F  /Vfc) 


(7-9) 


where  u  =  (uR,  u^)  ,  uR  is  the  rigid  body  control,  u^,  is  the  vector  of 

elastic  body  controls,  I  is  moment  of  inertia  of  the  undeformed 
structure,  is  the  vector  of  elastic  angular  momentum  coupling 
coefficients,  Mnn  is  the  modal  mass  matrix,  Knn  is  the  modal  stiff¬ 
ness  matrix,  F  is  the  modal  control  influence  matrix,  jr  =  [l,»»»,l]T 
(nxl ) , 


=  system  mass  matrix 


=  system  stiffness  matrix 


7-7 


j  =  system  control  influence  matrix 

U  p  J 

In  Eq.  (7-9),  0(tf)  and  n(tf)  are  prescribed  and  0(tf), 
j{(tf),  UR(tf),  ug(tf)  are  unknown.  Thus  to  obtain  a  solution 
for  u^(tf),  the  following  dynamical  constraints  are  imposed  on  0(tf) 
and  r[(  tf ) ,  name  ly : 


0 ( t  )  =  0  .(t_) 

f  target  f 


(7-10) 


n(tf)  =  o 


(7-11  ) 


As  a  result,  the  rigid  body  motion  osculates  the  target  motion  through 
two  time  derivatives,  and  the  flexible  body  response  is  zero  through  two 
time  derivatives,  where  Sf  is  given  by 


s,  ■=  [0  ( t .  )  0«»»0  0  0»*»0l. 

— f  L  target  f  taroet  f 


target  f ' 


Subject  to  the  given  boundary  conditions  of  js(tf)  and  the 
constraints  of  Eq.  (7-10)  and  (7-11),  the  solution  for  u(tf)  in  Eq. 
(7-9)  follows  as 


u(t£)  =  P  0„  fc(t.) 

—  f  target  f 


(7-12) 


where 


*  T  -1  T 

P  =  (B  WB)  B  W 


w  denotes  the  least-squares  weighting  penalty  on  the  rigid 

body  accelerations 

W  for  i  =  1,*«»,n  denotes  the  least-squares  weighting 

Vi 

penalty  on  the  i-th  modal  amplitude  accelerations 


1  4r 


=  control  influence  matrix 


In  Eq.  (7-12),  P*  is  the  weighted  least-squares  pseudo-inverse  operator 

corresponding  to  B"1 ,  since  B  is  in  general  rectangular. 

Assuming  that  the  spacecraft  control  system  is  turned  off  at  the 
start  of  the  retargeting  maneuver,  u^O)  is  given  by 


u  (0)  =  0 

—a  — 


(7-13 


Substituting  Eq.  (7-12)  and  (7-13)  into  Eq.  (7-8),  leads  to 


,  4>  Mt  ) 

s  X  sub  1  _  o 


$  .  $  II  ) 

V  VbJl*-*  ° 


s(t  )  -  <j)  s(t  ) 

—  f  'S3 —  O 


(7-14 


P*  ..  (t.)  -  6  s ( t  ) 

target  f  ua8  —  ° 


Equation  (7-14)  is  in  the  standard  form  for  the  linear  algebraic 
equation  Ax  =  Ja,  which  can  be  solved  via  Gaussian  elimination  for  jt, 
yielding  \(tQ)  and  Ub(tQ). 


The  optimal  control  time  histories  for  the  retargeting  maneuver 
are  obtained  by  integrating  Eq.  (7-5),  (7-6),  and  (7-7)  subject  to  the 
initial  conditions  given  by  s(t  ),  \(t  ),  and  U(t  )  =  [0T  ,  u?(t  )  1T. 

~  O  ”■  O  —  o  L— ■  — D  O  J 


7,4  Optimal  Integral  Pre-Compensation  LQ  Regulator  Design 

In  this  section  the  open-loop  control-rate  penalty  technique 
presented  in  Ref.  7-1  is  extended  to  the  closed-loop  case.  In 
particular,  for  the  results  of  this  section  it  is  assumed  that  the 
performance  index  penalizes  the  states,  controls,  and  first  time- 
derivatives  of  the  controls. 

The  first  step  in  the  problem  formulation  consists  of  adding  the 
control  to  the  state  equation  given  by 


As  +  Bu  ( 2Nx1  ) 


(7-15) 


where  £  is  the  state,  \i  is  the  control,  A  is  the  system  dynamics  matrix, 
and  B  is  the  control  influence  matrix.  As  a  result,  the  augmented  state 
equation  becomes 


•  tm  m 

Z  *  KZ  +  Bu2  [ ( 2N+NC ) xl ]  (7-16) 

where 
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0 

j5 

,  B  = 
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L  _ 

u,  -  u  (N  xl )  u_  =  u,  (N  xl ) 

— i  —  c  —2  — 1  c 


and  Nc  is  the  number  of  independent  controls 


7.4.1  Optimal  Control  Problem 


The  optimal  control  problem  is  to  seek  a  solution  of  Eq.  (7-16) 
which  minimizes 


J  =  J  >  9>t  I<tf )  +  j  ftf  [lTQl  +  £2 R  u2]dt  (7-17) 


subject  to  prescribed  initial  conditions 


£(t  ) 
—  o 


s(t  ) 
—  o 

u,  (t  ) 
— 1  o 


where  £Pf  is  the  augmented  terminal  state  weight  matrix,  Q  is  the 
augmented  state  weight,  and  R  is  the  weight  matrix  for  the  control- 
rate.  As  shown  in  what  follows,  the  commanded  control  becomes  the 
control  rate,  and  the  solution  follows  as  a  standard  LQ  regulator 
design. 

The  initial  control,  jii(t0),  in  Eq.  (7-17)  can  be  prescribed 
since  the  resulting  commanded  control  satisfies  a  first  order  differ¬ 
ential  equation  with  a  free  boundary  condition.  As  shown  in  Ref.  7-1 , 
it  is  convenient  to  select  _ui(t0)  *  since  this  choice  minimizes 
the  initial  control  jump  discontinuity  which  can  excite  the  residual 
plant  dynamics  of  flexible  space  structures. 

From  Pontryagin's  principle  it  easily  follows  that  the  necessary 
conditions  for  the  augmented  state,  augmented  costate,  and  control-rate 
are  given  by 


State: 


E 


AE  -  BR_1BTA 


(7-18) 


Costate: 


(7-19) 


Control-Rate:  u? 


-R"1BtA 


(7-20) 


7.4.2  The  Integral  Pre-Compensation  Regulator  Design 

Since  the  final  augmented  states  are  free,  Pontryagin's  principle 
provides  the  following  terminal  boundary  condition  relating  the 
augmented  states  and  costates 

A  (tf)  =  _£(tf)  (7-21) 


As  a  result,  it  can  be  shown  that  a  feedback  solution  can  be  obtained  by 
letting 


A(t)  =  ^(t)E(t) 


(7-22) 


for  t  e  [tQ,  tj.]. 

Upon  taking  the  time  derivative  of  Eq.  (7-22)  and  introducing 

•  • 

1  and  K  from  Eq.  (7-18)  and  (7-19),  one  finds  that  the  matrix  ^(t) 
satisfies  the  following  differential  matrix  Riccati  equation 

ff  =  -<jPa  -  -  Q  (7-23) 

which  is  subject  to  the  boundary  condition  given  by 


sn  tf)  = 


(7-24) 


where  is  defined  in  Eq.  (7-17). 

Substituting  Eq.  (7-22)  into  Eq.  (7-20)  yields  the  required 
control-rate  command  as  the  following  linear  combination  of  the 
augmented  states: 


7-12 


U2(t)  -  -«_1BT.^(t)£(t) 


(7-25) 


By  partitioning  the  matrix  ,r/(t)  as  follows: 


3>  = 


^12  I  }  2N 


^21  ^22  }  Nc 


<y  =  ^ 

21  12 


(7-26) 


and  carrying  out  the  implied  matrix  multiplication  in  Eq.  (7-25),  the 
control-rate  command  can  be  written  as: 


where 


u  (t)  =  -  Gs(t)  -  G  u,(t) 

—2  —  u — 1 


G  *  R  .^21(t) 


(7-27) 


G„  -  *  ■/'22{t) 


From  Eq.  (7-27)  it  follows  that  the  control-rate  command  is 
obtained  by  feeding  back  the  current  states  and  controls  through  time- 
varying  gains.  Alternatively,  if  the  infinite  time  problem  is  of 

• 

interest,  then  in  Eq.  (7-23)  is  set  to  zero  and  the  optimal  gains  are 
obtained  from  the  solution  for  the  algebraic  matrix  Riccati  equation: 

.^A  +  AT,9?  -  .^BR"1  BT,9f’  +  Q  =  0 


7-13 


The  block  diagram  for  the  integral  pre-compensation  LQ  regulator 
design  is  shown  in  Figure  7-2. 

From  Eq.  (7-27)  it  follows  that  since  U2  integrated  to 
provide  ^i(t),  one  can  interpret  the  control  scheme  of  this  section  to  be 
an  integral  pre-compensation  strategy.  However,  given  the  high- 
frequency  roll-off  characteristics  of  this  approach,  Eq.  (7-27)  can  also 
be  validly  interpreted  as  a  low-pass  filter  for  the  control.  Of  special 
interest  is  the  observation  that  a  single  solution  of  the  Riccati 
equation  yields  both  the  regulator  gains  and  the  low-pass  filter  gains. 

The  results  of  this  section  can  be  generalized  for  the  case  of  n-th 
order  time  derivatives  appearing  in  the  performance  index,  leading  to  a 
k-th  order  commanded  control  of  the  form 


u  =  -G  u-G»u-...-G  u , _  .  —  Gs 

u-  u-  “(k-, )  >  - 


.here  (*)(k)  S  ^  (•) ,  G  -  *'Vk+,  , 

at 

G  =  R-1^.  ,  .  ,  G«  =  R*1^.  ,  G 

u  k+1,2  u  k+1, 3  u 


(k-1 ) 


R_1  SP, 


k+1 ,k+1 


sp  = 


SP,A 

SP 

...  SP 

11 

12 

1 

SP 

SP 

...  SP 

21 

22 

2 

Cf  cp 

'k+1  , 1  *k+1,2 


2, k+1 


...  .9 

k+1 ,k+1 


The  block  diagram  for  the  k-th  order  integral  pre-compensation  LQ 
regulator  design  is  shown  in  Figure  7-3. 


After  some  reductions,  the  block  diagram  in  Figure  7-2  can  be 
written  as  shown  in  Figure  7-4.  The  loop  transfer  matrix  broken  at  the 
plant  input,  (a),  is  given  by 

“loop'"’  *  (SI>!  *  Gu,''G(8I2H  -  A)''B  <7-28 

c 

where  I  =  diag( 1 , • • • , 1 )  is  (N  xN  )  and  I_„  =  diagd , •••,1 )  is  (2Nx2N) 
N  C  C  2N 

C 

From  Eq.  (7-28)  it  follows  that  Gi00p(s)  has  the  following 
properties : 

(1)  A  two  pole  roll-off  in  each  channel; 

(2)  The  poles  of  G^00p(3)  consist  of  the  original  plant  poles, 
plus  the  eigenvalues  of  -Gu;  and 

(3)  For  scalar  input,  Gu,  as  determined  from  the  Riccati 
equation,  is  a  break  frequency. 

The  high  frequency  roll-off  characteristics  of  Gi00p ( s )  are 
believed  to  be  particularly  significant  for  the  residual  plant  dynamics 
of  flexible  space  structures. 


7.4.4  Observer-Based  Compensation  Design 


At  the  present  time  the  following  two  approaches  are  being 
considered  for  generating  the  observer  gains  K  required  in  the  state 
estimator.  The  first  approach  employs  a  standard  root  square  locus 
technique.  The  second  approach  employs  the  Doyle-Stein  loop  transfer 
recovery  procedure  to  improve  the  crossover  properties  of  the  loop 
transfer  function  broken  at  the  plant  input. 


The  Doyle-Stein  procedure  consists  of  solving  the  Kalman-Bucy 
filter  (KBF )  Riccati  equation 


where 


Mo(q> 


o  -  utn'tilql  -  ncTN_1cn 

O 


2  T 
M  +  q  B  B 
O 

nominal  plant  noise  intensity  matrix 


(7-29) 


=»  observation  noise  intensity  matrix 
=  Doyle-Stein  robustness  parameter 


and  the  observer  gains  are  given  by 


T  -1 

K  =  -  n  CN 


(7-30] 


As  q  ♦  ®,  in  the  KBF  Riccati  equation  above,  the  loop  transfer  matrices 
T^(s)  with  the  loop  broken  at  point  i  in  Figure  7-5,  have  the 
following  significant  limiting  properties 


T2(s)  ♦  T1  (s)  as  q  >  « 


(7-31 


where 


T1  (s)  =  $2(s)  G  $  (s)  B 
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The  robustness  recovery  procedure  is  attractive  because  T-j(s) 
has  the  standard  LQ  stability  robustness  guarantees: 


amin f1  +  —  1  for  a^1  “  — 


(7-32) 


and 


°min t 1  +  T i1  (ju>  ]  —  "2  for  a11  “ 


(7-33) 


where  o^n  denotes  the  minimum  singular  value.  In  particular,  T2(s) 
is  constrained  to  retain  a  4-pole  roll-off  at  high  frequency,  where  two 
poles  are  contributed  by  the  plant,  one  by  the  integrator,  and  one  by 
the  rest  of  the  compensator.  In  practice  the  design  process  consists  of 
a  trade-off  between  the  low-frequency  disturbance  rejection  properties 
and  location  of  the  break  frequency  for  the  high-frequency  4-pole  roll¬ 
off.  In  the  control  system  design,  the  weight  matrices  in  the 
performance  index  of  Eq.  (7-17)  and  robustness  recovery  parameter  q  are 
the  design  variables  employed  in  the  trade-off  study. 


7.4.5  Closed-Loop  System  Stability 

Since  truncated  dynamical  models  are  used  to  characterize 
flexible  spacecraft  and  generate  preliminary  control  designs,  it  is 
important  to  assess  the  stability  of  the  resulting  control  designs  in 
the  presence  of  unmodeled  system  dynamics.  To  this  end,  the  equations 
governing  the  system  state,  control,  and  state  estimation  are  summarized 
as  follows: 


n 
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The  algebraic  sign  of  the  real  part  of  the  eigenvalues  of 
provides  the  desired  system  stability  information.  In  particular,  if 
the  algebraic  sign  of  all  of  the  real  parts  of  the  eigenvalues  of  is 
negative  or  zero,  then  the  system  is  stable.  Alternatively,  if  any  of 
the  eigenvalues  of  have  positive  real  parts,  then  the  system  is 
unstable.  Clearly,  if  instability  is  indicated  in  the  eigenvalue 
analysis  of  ,  then  the  control  design  must  be  modified  so  as  to  over¬ 
come  the  difficulty. 


7.5  A  Closed  Form  Solution  for  the  Differential  Matrix  Riccati 
Equation  for  Fixed  Time  Slewing  Maneuvers 

The  matrix  Riccati  differential  equation 

P(t)  =  -  ATP(t)  -  P(t)A  +  P(t)  BR_1BTP(t)  -  Q  (n  x  n)  (7-42) 

arises  in  the  solution  for  the  optimal  slewing  control  problem  for  the 
linear  time-invariant  system 

x(t)  *  A  x(t)  +  Bu(t),  x(t  )  *  x  (n  x  1 )  (7-43) 

—  —  —  o  — o 
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where  the  performance  index  is  given  by 


J  “  i  xT(t  )Sx(t  )  +1  /  f  [xT(t)Qx(t)  +  uTRu]dt  (7-44) 

2  £  1  2  o 

and  x(t)  is  the  state,  jj(t)  is  the  control,  A  is  the  system  dynamics 
matrix,  B  is  the  control  influence  matrix,  S  is  the  terminal  state 
weight  matrix,  Q  is  the  state  weight  matrix,  and  R  is  the  control  weight 
matrix.  If  R,  Q,  and  S  are  selected  to  be  symmetric  and  positive 
definite,  then  P(t)  is  symmetric  and  positive  definite.  The  positive 
definiteness  of  P(t)  is  important  for  the  results  that  follow. 

The  optimal  control  of  Bq.  (7-43)  is  the  following  linear  state 
feedback  law 


u(t)  -  -  R^B^miiU)  (7-45) 

where  P(t)  is  the  solution  for  Bq.  (7-42)  with  the  terminal  boundary 
condition 


P(tf)  -  S  (7-46) 

This  section  presents  a  simple  change  of  variables  for  B).  (7-42) 
which  permits  a  closed  form  solution  to  be  obtained  for  P(t),  when  R,  S, 
and  Q  are  selected  to  yield  P(t)  positive  definite. 

I 


p(t) 


(7-47) 


p  +  z-1  (t) 

ss 


where  Paa  is  the  solution  to  the  algebraic  matrix  Riccati  equation 


-A  P  -PA  +  PEP  -  Q  =  0  , 

ss  ss  ss  ss 


(7-48) 


Z(t)  is  a  matrix  function  to  be  determined,  and  E  =  B  R-1BT. 

The  substitution  in  Eq,  (7-47)  has  been  suggested  by  Potter  and 
VanderVelde  in  Ref.  7-2,  and  is  a  natural  consequence  of  the 
corresponding  results  for  scalar  Riccati  equations.  The  results  of  this 
section  extend  the  work  of  Potter  and  VanderVelde  by  providing  an 
efficient  analytic  solution  for  Z(t).  The  complete  solution  for  Z(t) 
can  be  found  in  Ref.  7-3. 

After  some  algebra  the  linear  constant-coefficient  matrix 
differential  equation  for  Z(t)  follows  as: 


Z(t)  =  G  Z(t)  +  Z(t)  G  -  E 


(7-49) 


where  G  =  A-EPgg,  and  G  is  the  steady-state  closed-loop  system  matrix 
possessing  eigenvalues  with  all  negative  real  parts. 

The  solution  for  Z(t)  is  obtained  by  introducing  a  coordinate 
transformation  which  diagonalizes  the  homogeneous  part  of  Z.  To  this 
end,  G  and  GT  in  Eq.  (7-49)  are  expressed  as 


G  =  R  A  L 


(7-50) 


T  T 

G  -  L  A  R 


(7-51 ) 


L 


where  R  is  the  matrix  of  right  eigenvectors  of  G,  L  is  the  matrix  of 
left  eigenvectors  of  GT,  A  is  the  diagonal  matrix  of  eigenvalues  of  G, 
and  R  and  L  are  normalized  as  follows 


T  T  — 1  T  — 1 

R  L  =  I  (R  =  L  and  L  =  R  ) 


Substituting  Eq.  (7-50)  and  (7-51)  into  Eq.  (7-49)  and  pre¬ 
multiplying  the  result  by  LT  and  post-multiplying  by  L,  leads  to 


T(t)  =  A  T(t)  +  T(t )  A  -  E 


(7-53 


where 


T(t)  «  Lt  Z(t)  L  (Z(t)=  R  T(t )  Rt) 
E  =  Lt  E  L 


Since  the  homogeneous  part  of  T{t)  is  diagonal,  the  solution  for 
T(t)  is  easily  obtained,  and  after  some  algebra  it  can  be  shown  to  be 


A(t~t  )  A(t-t  ) 

T(t)  =  e  [T(tQ)  “  D(t,tQ)]e  ° 


where 


=  diag  [A1  •••  An] 


A( t-t  )  A, (t-t  )  A  (t-t  ) 

i  =  diag  [e  •••  e  J 


-  P, .(t-t  ) 


initial  condition  matrix  for  T(t) 


T(t  )  = 

o 


[*]„  =  denotes  i-j  element  of  [*] 

As  a  result,  the  solution  for  Z(t)  is  known  and  solution  for  P(t)  can  be 
written  as 


P(t)  =  p  +  [R  T(t)  rT]_1  (7-55) 

ss  J 

The  initial  condition  matrix  for  T(t0)  is  evaluated  by  imposing  the 
boundary  condition  that  P(tf)  =  S.  After  carrying  out  the  necessary 
algebra  the  solution  for  T(t)  follows  as: 


[T<t)lij 


e 


where  H  *  LT(S  -  P  )  V 

ss 


(7-56) 


Using  Eq.  (7-55),  the  integral  pre-compensation  control  of  Section 
7.4.2  can  be  solved  for  fixed-time  slewing  maneuvers. 
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Figure  7-1.  Geometry  for  slewing  to  engage  and  track  a 
moving  target. 
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Figure  7-2.  Block  diagram  for  the  integral  precompensation 
LQ  regulator  design  using  full-state  and  full- 
control  feedback. 
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Figure  7-4.  Block  reduction  for  the  integral  precompensation 
LQ  regulator  design. 
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